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ABSTRACT 

We have combined data from the NOAO Deep Wide-Field Survey in Bootes and the Spitzer Space 
Telescope to determine basic properties for sixteen optically "invisible" MIPS 24 /xm (OIMS) and 
thirty-five optically "invisible" radio (OIRS) sources, including their spectral energy distributions 
(SED) and luminosities. Most OIMSs possess steep power-law SEDs over A res t = 1 — 10 (im, in- 
dicating the presence of obscured AGN in agreement with Spitzer spectroscopy. These objects are 
extremely luminous at rest-frame near and mid-IR (vL v (5 fim) ~ 10 38 — 10 39 W), consistent with 
accretion near the Eddington limit and further implying that they host buried QSOs. The majority 
of the IRAC detected OIRSs have flat 3.6 to 24 /mi SEDs, implying comparable emission from stellar 
photospheres and hot AGN illuminated dust. This may reflect relatively small amounts of dust close 
to the central engine or current low mass accretion rates. A small subset of OIRSs appear to be star- 
burst dominated with photometric redshifts from 1.0 to 4.5. The OIMSs and OIRSs with significant 
starburst components have similar Lk and stellar masses (M* as 10 11 Mq) assuming minimal AGN 
contribution. Roughly half of the OIRSs are not detected by Spitzer 's IRAC or MIPS. These are most 
likely z > 2 radio galaxies. The IRAC detected OIRSs are more likely than OIMSs to appear non 
point-like in the 3.6 \im and 4.5 fim images, suggesting that interactions play a role in triggering their 
activity. The AGN powered OIMSs may represent sub-millimeter galaxies making the transition from 
starburst to accretion dominance in their evolution to current epoch massive ellipticals. 
Subject headings: galaxies: high-redshift - galaxies: starburst - galaxies: active - infrared: galaxies 
- radio continuum: galaxies 
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1. INTRODUCTION 

Sensitive multi-band surveys in the optical and near- 
infrared portion of the electromagnetic spectrum have 
proven extremely successful in identifying large and di- 
verse galaxy populations in the high-redshift universe. 
Deep images obtained through optical U n GR filters, for 
example, readily find z > 3 galaxies by the presence of 
a strong 912 A Lyman break (e.g., Steidel et al. 1996). 
Lyman Break Galaxies (LBGs) tend to be massive galax- 
ies with star formation rates (SFR) between w 4 — 25 
Mq yr _1 and relatively low dust obscuration. Similarly, 
unobscured star forming disk (S-BzK galaxies) and pas- 
sively evolving spheroidal galaxies at z ~ 1.4 — 3 can be 
discerned from the myriad background sources by their 
optical/near-infrared colors or other magnitude-color cri- 
teria (see, e.g., Daddi et al. 2004; Adelberger et al. 
2004). The S-BzK galaxies in particular tend to possess 
large stellar masses (M» w 10 11 Mq) with SFR w 200 
M© yr _1 . Nevertheless, it had been long suspected that 
a substantial fraction of luminous z > 1 galaxies pow- 
ered by star formation or accretion onto super-massive 
black holes would not be identified using these and re- 
lated techniques because of heavy dust obscuration, but 
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would turn up instead at infrared, sub-millimeter, and 
radio wavelengths (see, e.g., Chary & Elbaz 2001 and ref- 
erences within). Indeed, deep but small-area surveys at 
sub-millimeter wavelengths by SCUBA (e.g., Smail et al. 
1997) revealed a population of extremely dusty galaxies 
at z ~ 1 — 3, whose high luminosities are generally though 
to originate primarily in star formation (SFR > 1000 M 
yr _1 ), though a small fraction may be AGN dominated. 
Galaxies powered by starbursts or active nuclei (AGN) 
can also be strong emitters at radio wavelengths through 
either large populations of supernova remnants or ac- 
cretion disk phenomena. Even though this emission is 
unaffected by obscuration, roughly 10 — 15% of compact 
radio sources identified in deep radio surveys have either 
extremely faint optical counterparts or none at all (e.g., 
Richards et al. 1999; Fomalont et al. 2002), and may 
represent distant luminous obscured starburst or AGN 
powered galaxies. Such objects make a still unknown 
contribution to the luminosity, chemical, and accretion 
history of the universe. Determining an accurate census 
of heavily obscured sources - their redshifts, luminosi- 
ties, space densities, and dominant power source - is one 
of the key issues in observational astrophysics and a pri- 
mary objective of the Spitzer Space Telescope (Spitzer, 
Werner et al. 2004). 1 

A number of ambitious large area multi-wavelength 
surveys incorporating Spitzer data have been conducted 
to date. These include the Great Observatories Origins 
Deep Survey (GOODS), the Spitzer Wide-area InfraRed 
Extra- galactic survey (SWIRE), and the First Look Sur- 

1 The Spitzer Space Telescope is operated by JPL, Califor- 
nia Institute of Technology for the National Aeronautics and 
Space Administration. Information on Spitzer can be found at 
http: / / ssc. spitzer . caltech. edu / . 



2 



Higdon et al. 



vey (FLS). This paper is concerned with the nature of 
infrared and radio selected sources with extremely faint 
optical counterparts originally selected from the NO AO 
Deep Wide-Field Survey (NDWFS) in Bootes (Jannuzi 
& Dey 1999) 2 with the expectation that they are either 
heavily extincted, at high redshift, or both. Higdon et 
al. (2005; hereafter, Higdon05) identified thirty-six Opti- 
cally "Invisible" Radio Sources (hereafter, OIRSs) out of 
377 compact or unresolved radio sources found in a VLA 
A-array 20 cm survey covering 0.5 deg 2 in the NDWFS 
Bootes field. These observations reached a flux density 
limit of « 80 /iJy (5a) at the three overlapping pointing 
centers. To be considered optically "invisible", a radio 
source must have no Byy, R, or / counterpart within a 
1.5" radius. This corresponds to limiting 3a magnitudes 
of approximately 26.9, 25.6, and 24.9 (Vega; 2" aperture) 
respectively, with the precise value depending upon loca- 
tion within the optical survey. The NDWFS region was 
surveyed at 24 /im and 70 /im with the Multiband Imag- 
ing Photometer for Spitzer (MIPS, Ricke et al. 2004). 
Analysis of the OIRS's (^parameter led Higdon05 to con- 
clude that they are a population powered by relatively 
unobscured radio-loud active nuclei. 3 While none of the 
OIRSs have measured redshifts, Higdon05 argued that 
they are likely to be at z > 1 based on the faint op- 
tical limits set by the NDWFS survey. An indepen- 
dently defined sample of seven optically faint (R > 24.5) 
and ten optically "invisible" sources from the NDWFS 
Bootes region, subject to the additional constraints that 
Fnftm > 0.75 mJy and vf„(24fim)/vf v (0.8fj,m) > 100, 
was observed using Spitzer 's Infrared Spectrometer (IRS, 
Houck et al. 2004). These Optically "Invisible" MIPS 
Sources (hereafter, OIMSs) are a high-z population, with 
z ~ 1.6 — 2.7. Comparisons with mid-IR spectra of local 
starburst and AGN dominated galaxies suggested that 
the primary energy source in thirteen of the seventeen 
OIMSs is a heavily obscured active nucleus (Houck et al. 
2005, hereafter Houck05). All seventeen have inferred 
8 — 1000 /im luminosities of wlO 13 L Q , placing them in 
the "hyper" -luminous class. 

It is not known how the radio and infrared selected 
"invisible" populations in Higdon05 and Houck05 are 
related, apart from the fact that the majority of both 
appear to possess powerful AGN. We do not know, for 
example, if the OIRS population lies at systematically 
higher redshifts than the OIMSs with similar Lir, or for 
that matter, if they are a sub-L* population at z ~ 0.5. 
Key to addressing these issues is estimating photometric 
redshifts for the OIRSs, which are generally too faint for 
Spitzer's IRS. In this paper, we present results obtained 
by combining optical B\y, R, and I, and MIPS 24 /im, 
70 /im, and 160 /im flux densities (or limits) with data 
from Spitzer's Infrared Array Camera (IRAC, Fazio et al. 
2004) Shallow Survey of the NDWFS Bootes field (Eisen- 
hardt et al. 2004). The IRAC Shallow Survey covers 8.5 

2 The NOAO Deep Wide-Field Survey is supported by the 
National Optical Astronomy Observatory, which is operated by 
AURA, Inc., under a cooperative agreement with the National Sci- 
ence Foundation. Information on the NDWFS can be found at 
htt p : / / www .noao.edu/noao/noaodeep/. 

3 q = /og(F24 Mm /-F20cm)- Systems powered by star formation 
or radio-quiet active nuclei possess q = 0.5 — 1.1. Smaller, and in 
particular, negative values of q indicate the presence of increasingly 
radio-loud active nuclei (Appleton et al. 2004). 



deg 2 in four bands centered at 3.6, 4.5, 5.8, & 8.0 /zm. 
For z w 1 — 3, these bands measure emission in the rest- 
frame near-infrared. Because this wavelength regime is 
much less sensitive to extinction relative to the optical, 
the IRAC data can constrain the evolved stellar content 
of galaxies at this epoch. Photometric redshifts are also 
possible by virtue of the rest-frame 1.6 /im peak arising 
from the H opacity minimum in the photospheres of 
evolved stars. It is our aim to derive basic properties for 
both optically "invisible" source populations, including 
rest- frame near- infrared luminosities and SEDs, and to 
determine, for example, if their rest-frame near-infrared 
emission is dominated by starlight or an AGN. Only then 
would it be possible to relate OIMSs and OIRSs to other 
high-z populations. It is also our intent to use the IRAC 
images to explore the near environments of OIMSs and 
OIRSs for additional clues to their nature and evolution. 

A brief description of the extraction of IRAC flux 
densities for the infrared and radio selected "invisible" 
sources is given in §2. Photometric redshift estimates 
for selected OIRSs are presented in §3, along with deter- 
minations of the average rest-frame SEDs, luminosities, 
and near environments of the OIMS and OIRS samples 
detected with IRAC. The nature of these two popula- 
tions is discussed in §4. These results are summarized in 
§5. Throughout this paper we will refer to the OIRSs and 
OIMSs by number (i.e., OIRS #97 or OIMS #13), corre- 
sponding to their entries in Table 2 of Higdon05 and Ta- 
ble 1 of Houck05. Source coordinates can be found from 
those papers directly, or through their SIMBAD desig- 
nations, [HHW2005] # and [HSW2005] # for the OIRSs 
and OIMSs respectively. We assume a flat ACDM cos- 
mology, with fiM=0.27, Cl\=0.73, and a Hubble constant 
of 71 km s" 1 Mpc" 1 . 

2. OBSERVATIONS AND SOURCE EXTRACTION 

Flux densities or limits at 3.6, 4.5, 5.8, and 8.0 /im 
for the OIRSs and OIMSs were extracted from the fi- 
nal IRAC Shallow Survey mosaics (pipeline version Sll) 
using SExtractor (Bcrtin & Arnouts 1996) with 4.8" di- 
ameter apertures. All sources are nearly point-like, and 
point-source aperture corrections were applied. The 3.6 
/jm images were used for matching the radio and 24 /im 
selected source positions since they had the highest point 
source sensitivity (F 3 ,a^ m — 6.4 /iJy (5a) in a 4" diam- 
eter aperture) and angular resolution (1.7" FWHM). To 
qualify as a "match" the coordinates of an OIMS or OIRS 
had to agree with an IRAC 3.6 /im source to within 2". 
We decided upon this match radius after comparing the 
radio and i*24^m > 0.75 mJy MIPS catalogs with the 
IRAC Shallow Survey catalog. For example, 2/3 of the 
compact 20 cm radio sources (excluding the OIRSs) had 
a 3.6 /im source counterpart centered within 0.65". Since 
the great majority of these matches are real we took 0.65" 
to represent the la relative positional accuracy between 
the radio and IRAC 3.6 /im catalogs. A similar value 
(0.71") was found for the la relative positional accuracy 
between the F24^ m > 0.75 mJy MIPS sources and 3.6 
/xm counterparts in the Shallow Survey. The match ra- 
dius of 2", while seemingly small, represents essentially a 
3a criterion, i.e., we would expect « 1 radio/IRAC mis- 
match out of the full 377 compact radio source sample. 
As a consequence, OIMSs and OIRSs can be said to lack 
IRAC counterparts at the 3a level if no 3.6 /im sources 
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are within 2", and at the 5a level if there are none within 
3.3". 

The original count of OIRSs in Higdon05 was thirty- 
six. However, we are no longer confident that #362 is a 
clear-cut OIRS. As a result, we will not include it in the 
following analysis. Of the remaining thirty-five OIRSs 
within the IRAC Shallow Survey, nineteen (54%) were 
unambiguously detected at 3.6 /im as shown in Figure 1. 
Even when a potential optical counterpart exists within 
~ 2 — 3" (e.g., #97 & 410) the compact radio source sits 
squarely on an IRAC source. Of these, six (17%) were 
detected in three IRAC bands to > 3a (OIRS #79, 176, 
208, 245, 346, and 349), and two (#176 and 349) in all 
four. The median 3.6 /an and 4.5 /zm flux densities for 
these nineteen objects are 18.3 /iJy and 22.5 /iJy, respec- 
tively. Sixteen OIRSs were found to lack IRAC counter- 
parts within 2" and are shown in Figure 2. Note that 
in four instances (#9, 305, 430, and 441) possible IRAC 
counterparts are present roughly 3 — 5" away. OIRSs #9, 
305, and 441 are particularly suspicious as their candi- 
date IRAC matches are also optically "invisible", i.e., 
like the 19 OIRSs with clear IRAC counterparts. The 
coordinate mismatches in these instances however equal 
or exceed the radio/IRAC relative positional accuracy at 
the 5a level. Moreover, the expected number of OIRSs 
in Figure 2 with an optically "invisible" IRAC source 
3 — 5" away due to chance is w 3, given the estimated 
density of 3.6 /im sources with no optical counterparts 
(3.3 arcmin~ 2 ) in the Shallow Survey. We conclude that 
these four OIRSs are unlikely to be physically associated 
with the relatively nearby IRAC sources. However, even 
if we allow these four matches, it remains true that a 
significant fraction (at least 12/35 or 34%) of the OIRSs 
lack clear 3.6/zm counterparts in the IRAC Shallow Sur- 
vey. 

It should also be noted that while the radio source 
morphologies in Figures 1 and 2 are for the most part 
quite simple, extended radio emission may have been 
lost due to a lack of short spacings (though the 4.5-hours 
spent at each pointing center in the A-array observations 
should reliably represent emission up to 12" scale-sizes) 
or (1 + z)~ A surface brightness dimming that might oth- 
erwise have connected an OIRS to an optical or infrared 
source. Observations in more compact array configura- 
tions will be needed to resolve this issue. With these 
caveats in mind, stack averages of the 3.6, 4.5, 5.8, and 
8.0 /im images at the radio source positions led to flux 
density upper-limits of 1.5, 2.0, 10.1, and 11.0 /iJy (3cr) 
for the sixteen OIRSs with no apparent IRAC counter- 
parts. 

OIMS #11 was not included in this study due to its 
uncertain redshift. 4 All of the sixteen remaining OIMSs 
were detected by IRAC, although OIMSs #17 was only 
observed at 3.6 /zm and 5.8 /zm. These are shown in 
Figure 3 and their flux densities or limits are listed in 
Table 2. The median 3.6 /im and 4.5 /im flux densities 
for the fifteen OIMSs observed in all four IRAC bands 
are 17.0 /zJy and 26.1 /zJy, i.e., quite similar to that of the 
nineteen IRAC detected OIRSs. However, a much higher 
fraction of these OIMSs (12/15 or 80%) were detected in 

4 The S/N for this source's IRS spectrum is low. While Houck05 
gives Zapec = 0.70 ±0.24, a redshift of fs 2.5 is just as feasible using 
the same galaxy template. 



all four IRAC bands. The OIMSs are also significantly 
brighter at 8.0 /im than the OIRSs: the median OIMS 
Fg.ojmi is 136 /iJy, while the median i^.o^m for the nine 
OIRSs detected in this band is 37 /zJy. 

The OIMSs were chosen from the MIPS 24 /zm catalog 
to satisfy > 0.75 mJy, and the sources in Table 2 

have a median -F^m of 1.1 mJy. By contrast, only four 
OIRSs were detected at 24 /zm at 5a or higher (OIRS 
#97, 176, 245, and 363), with a mean F 2 ^ m of 0.38 
mJy. Two additional OIRSs were detected at the 3a level 
(OIRS #208 and 349), with a mean F 2 4 Mm of 0.24 mJy. 
We stack-averaged MIPS 24 /im sub-images centered on 
the eleven non-detected OIRSs in an attempt to recover 
faint emission from these sources. No significant signal 
was found, allowing us to set a flux density upper-limit 
of 67 /zJy (3a) at 24 /im for these sources as a group. No 
OIRSs or OIMSs were detected with MIPS at 70 /im or 
160 /im, placing 3a upper-limits of 24 mJy and 60 mJy 
at these two wavelengths. 

The IRAC and MIPS flux densities (or 3a limits) for 
the two source populations are listed in Tables 1 and 2. 
To these have been added flux densities (or limits) in 
B w , R and I from the NDWFS catalog 5 , and at 20 cm 
from cither Higdon05 for the OIRSs or de Vries et al. 
(2002) for the OIMSs. 

3. PROPERTIES OF THE OPTICALLY INVISIBLE 
POPULATIONS 

3.1. Photometric Redshifts 

Photometric redshifts (z p hot) were calculated using a 
set of galaxy templates corresponding to starburst (M 82, 
from Xu et al. 2001), composite starburst plus AGN 
(Arp 220, from Devriendt et al. 1999, modified using IRS 
spectra to better represent its polycyclic aromatic hydro- 
carbon (PAH) emission and silicate absorption), AGN 
dominated (I Zw 1, Xu et al. 2001), and elliptical (Silva 
et al. 1998) galaxies. For each OIRS, the templates were 
rcdshifted from z = 0.1 — 6.0 in steps of 0.05, re-gridded 
and scaled to the observed flux densities. Redshifts were 
determined by minimizing the reduced x 2 , defined as 
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(O a - a z T (1+z)A ) 2 



(1) 



where Oa and T( 1+z )a are the observed and template flux 
densities, Ndof are the degrees of freedom, and a z is a 
scaling constant calculated at each value of z. The opti- 
cal and far-infrared upper-limits were included in the \\ 
calculation only when they were exceeded by the tem- 
plate flux density (i.e., when a z T(i +z ^\ > Oa). When 
this occurred A was set to the respective limit and \\ 
was calculated as above. Formal uncertainties in z p hot 
were estimated by numerically evaluating 



2.0 



d 2 xl/dz 



(2) 



Additional extinction corresponding to Ay = 1-2 (assum- 
ing a screen geometry and the Galactic reddening law in 

5 If not already published, new 3cr limits were calculated using 
m — 2.5 log(3 o" p ; x ^/N p ; x ), where <x p i x is the rms per pixel, N p i x 
is the number of pixels within the 4" diameter aperture, and m a 
is the appropriate zero-point constant. 
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Mathis 1990) was typically required to satisfy the opti- 
cal constraints. 6 This was applied to the template in 
Equation 1 prior to redshifting. 

As a test, we determined photometric redshifts for the 
four OIMSs that Houck05 fit with cither an Arp 220 or 
NGC 7714 spectral template (OIMS #2, 7, 12, and 14 
in Table 2). The minimum-^ 2 fits are shown in Figure 4 
along with the corresponding photometric redshifts. We 
find good agreement between the photometric and spec- 
troscopic redshifts of all four, with |z p hot — z spcc \ = 0.1, 
0.1, 0.1 and 0.3, respectively. The Arp 220 template 
with an additional w 1.5 magnitudes of visual extinc- 
tion provided the best fits for all four, including OIMS 
#2, whose IRS spectrum was classified as a starburst 
(NGC 7714) by Houck05. Starburst templates gave sim- 
ilar Zphot for this source, but produced too little emission 
at 8.0 /an and 24 /im, implying the presence of an ad- 
ditional dust component in the rest-frame mid-IR. The 
fact that OIMSs #12 & 14 are significantly brighter than 
the Arp 220 template in the 24 /im band may reflect 
stronger PAH emission in these sources. Again, similar 
photometric redshifts were obtained using the M 82 tem- 
plate and down-weighting the 24 /im measurements, how- 
ever only the Arp 220 template could match the observed 
-p24/xm/-f8.0Atm- The remaining OIMSs in Table 2 either 
lacked a sufficient number of data points for template fit- 
ting (e.g., #5, 10, and 17) or possessed power-law SEDs, 
for which unambiguous redshift determinations were not 
possible. 

Photometric redshifts were derived for six OIRSs: #97, 
176, 208, 245, 349, and 363. Each was detected in at least 
three IRAC and MIPS bands to > 3er. The minimum-x 2 
fits are shown in Figure 5, along with the best match- 
ing template and derived z p hot- The results are summa- 
rized in Table 3. In each case we are fitting the inflec- 
tion point in the spectrum where the dominant emission 
source changes from evolved stars to warm dust. For 
OIRS #363 the inflection point is implied by the upper- 
limits in the IRAC 5.8 /im and 8.0 /im bands, which effec- 
tively rule out a power-law between 3.6 /im and 24 /im. 
The derived redshifts range from 1.0 to 4.5 with a mean 
of 2.2. The formal redshift uncertainty is w 0.2 — 0.6, 
though a comparable contribution to the error budget 
may arise from our ignorance of the source's true SED. 
For example, OIRSs #97 and 245 would have been de- 
tected at 160 /im to > 3a had the Arp 220 template 
accurately represented their far-infrared SEDs. (A more 
in depth analysis employing a larger set of templates and 
extinction laws is not warranted at this point given the 
paucity of data points.) The six OIRSs are best repre- 
sented by either the starburst M 82 (#176, 208, and 349) 
or the starburst/AGN composite Arp 220 (#97, 245, and 
363) templates. In each case the best fitting template 
was unambiguous. In particular, the AGN-dominatcd 
galaxy templates make poor matches to the SEDs in Fig- 
ure 5. Likewise, elliptical galaxy templates are ruled out 
by the detection of emission at 24 /im. This sub-sample 

6 Estimates of the intrinsic extinction in the template galaxies 
can vary considerably depending on angular resolution, precise lo- 
cation in the source, and method used. Typical values of arc 
80 magnitudes for Arp 220 (Devriendt ct al. 1999), 10 for the cen- 
tral region of I Zw 1 (Eckart et al. 1994), to 5 for M 82 (Lester et 
al. 1990; Telesco et al. 1991). The elliptical models of Silva et al. 
have essentially zero extinction. 



of OIRSs detected both with IRAC and with MIPS at 24 
/im represents a population of sources with a significant 
starburst component with a redshift range similar to the 
OIMSs. 

Two other OIRSs, #79 and 346, permitted rough 
photo-z estimates. Minimum-x 2 fits are shown in Fig- 
ure 6. For #79, the fiat IRAC spectrum could reason- 
ably be interpreted as the rest-frame ~l-3 /im region, 
consistent with redshifts between 1.5 — 2.5. In fact, the 
smallest x 2 is achieved with a starburst (M 82) template 
at z p hot~ 1-6 with one magnitude of additional extinc- 
tion. An Arp 220-like SED appears to be ruled out unless 
the 9.7 /im silicate absorption feature in these sources 
is much deeper than the template. We emphasize that 
these redshifts are approximate, as z p hot = 2.5 also yields 
a small x 2 - The second OIRS in Figure 6, #346, shows 
a fairly steep power-law shape (a = -2.1, for F v oc v a ) in 
the IRAC data that does not extend to the MIPS 24 /im 
band. If this spectral break occurs at rest-frame wave- 
lengths of 1-2 /xm, then redshifts exceeding 3 are implied, 
and a reasonable fit using the M 82 and Arp 220 template 
at Zphot = 3.5 — 4.6 are possible, as is shown in Figure 
6. Note that good fits also result using the elliptical 
template for both, since the 24 /tin upper-limit provides 
no constraint to the shape of the SED. In both cases, 
the derived z p hot is close to that derived using M 82 or 
Arp 220. Because these redshifts are only estimates, we 
do not include them in Table 3. 

Overall, no starburst, AGN, or composite template 
could be found that would permit a fitted z p hot of less 
than one for the above eight OIRSs. At the very least, a 
significant fraction of OIRSs are at high redshift. While 
galaxy templates with a sizable young stellar component 
(i.e., Arp 220 or M 82) are consistent with the data, 
early-type galaxies cannot be ruled out in the two cases 
where there are no 24 /im detections (Figure 6). 

There are insufficient data to determine photometric 
redshifts for the remaining OIRSs in Table 1. It is im- 
portant to keep in mind that we have derived reasonably 
confident Zphot for only 1/3 of the IRAC detected OIRSs. 
This represents only 17% of the total OIRS sample from 
Higdon05. We will address the issue of how likely these 
six are to be representative of the entire OIRS population 
below. 

3.2. Spectral Energy Distributions and Emission 
Mechanisms 

HouckOS found that the mid-IR spectra of the OIMSs 
could be fit using only four templates, all derived from 
local galaxies (see Table 4). The same templates also fit 
the wider spectral energy distributions examined here. 
They are: (1) a heavily obscured Seyfert 2 nucleus, repre- 
sented by IRAS F00183-7111 and denoted obsSy2, (2) an 
obscured Seyfert 1 nucleus, represented by Mrk 231 and 
denoted obsSyl, (3) a starburst nucleus, represented by 
NGC 7714 and denoted sb, and (4) a combination of an 
obscured AGN and a starburst, represented by Arp 220 
and denoted comb. The first four panels of Figure 7 show 
average rest-frame SEDs for the sources in each category. 
The SEDs were found by shifting each source's flux den- 
sities (or limits) to the rest-frame, then normalizing at 
A rcst = 1.4 /tm. The individual scaled SEDs were then 
coadded and binned. 

The averaged SEDs for the obsSyl and obsSy2 OIMSs 
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both show steep power-laws extending from A ros t ~ 1 — 
10 fim. This is consistent with a luminous and heavily 
obscured active nucleus. Assuming this emission arises 
from dust, which is reasonable given the high extinction 
deduced by Houck05, emission at ~1 /im implies that 
some of the dust approaches the sublimation temperature 
of refractory dust. For the obsSy2 SED, the power-law 
exponent a « —4.0 (F u oc z/*) over the range 1 < A rost < 
4 /<m, and there is a break in the power-law at A res t ~ 5- 
8 /im. For the obsSyl OIMSs, the slope is shallower 
(a rts —2.3), but the power-law component extends over a 
wider range, 1 < A rcs t ^ 10 fim, with no break. The slope 
is still steep compared to SEDs of unobscurcd AGN (e.g., 
Elvis et al. 1994) and even compared to the lcnscd z = 
3.91 quasar APM 08279+5255 (Soifer et al. 2004), for 
which a = —2.0 over 1 < A rcst < 4 /im. The steep SEDs 
suggest a distribution of dust mass that rises rapidly with 
distance from the central heating source, and the break 
in the obsSy2 SED suggests that the rapid rise ends at 
a distance corresponding to Tdust 400 K. The absence 
of PAH emission in the IRS spectra of these objects is 
consistent with the AGN interpretation of the emission. 

The sb SED rises about as steeply (a w —2.3) as the 
obsSyl SED in the rest-frame 1-8 [im spectral range, 
though with only one object the slope is poorly deter- 
mined. A 1.6 /im bump can be seen in the IRAC data, 
indicating that the rest-frame 1-3 /jm emission arises pri- 
marily from evolved stars. The steep rise at longer wave- 
lengths is characteristic of dust in star formation regions. 
In contrast, the average SED for comb objects is much 
shallower for A ros t between 1-3 fim, a w —1.1, with a 
steeper rise from rest 3 to 8 /im. The lack of a prominent 
1.6 /im bump suggests a power- law AGN component con- 
tributes much of the rest-frame 1 — 3 /im flux, and the 
steepening suggests the presence of dust associated with 
the star burst, consistent with the dual starburst/AGN 
nature of these objects. 

The average SED for the six IRAC-detected OIRSs 
with photometric redshifts is also shown in Figure 7. 
This SED resembles the sb class by virtue of the clear 
1.6 /im bump. However, the steep rise for A ros t between 
3-10 /im may reflect a combination of star formation 
plus a weak active nucleus. In any event, there is no 
evidence of a dominant obscured active nucleus, which 
suggests that a significant fraction of the luminosity is 
provided by star formation. This is in marked contrast 
with the OIMSs, which are powered in most cases by 
dusty AGN. We again wish to stress that the averaged 
OIRS SED was created using only the six sources with 
photometric redshifts and may not be representative of 
the remaining thirteen IRAC detected OIRSs in Table 1, 
much less the parent population of OIRSs in Higdon05. 

Additional insight into the dominant power source for 
OIMSs and IRAC detected OIRSs can be gathered from 
their positions in an IRAC two-color diagram. Figure 
8 shows that OIMSs and OIRSs tend to occupy dis- 
tinct regions. Most OIMSs define a fairly narrow dis- 
tribution extending from relatively flat observer-frame 

SEDs (i'8.0 J um/-F4.5 J um = F 5 ,Sfim/ F^Sfim ~ 1) in the 

lower-left to steep power-law SEDs (Fg.o/im/^i.s^m = 
F 5 .8nm/F 3 . 6f im ~ 10) in the upper-right. Nearly all of the 
obsSyl and obsSy2 OIMSs lie above fs.o^m/^k.Sjum = 3, 
in the portion of the diagram dominated by quasars and 
obscured AGN (Lacy et al. 2004). This is consistent 



with both their IRS spectra (which show deep 9.7 /im 
silicate absorption and no PAHs) and their steep power- 
law SEDs in Figure 7. The comb and sb OIMSs are situ- 
ated at F 8 .o^m/F 4 . 5lim < 3 and F 5 . 8Aim /F 3 . 6Aim < 2, where 
starburst systems dominate, also in agreement with their 
IRS spectra. 

The IRAC detected OIRSs are shown in the lower 
panel of Figure 8. The six OIRSs with photometric 
redshifts have positions that are for the most part con- 
sistent with the dusty M 82 track for z w 1 — 3, in 
agreement with their z p hot and fitted spectral templates. 
The remainder of this population are also clustered in 
the lower-left of the figure, i.e., -F8.o J um/-F4.5 jUm < 3 and 
F5.8fj.rn/ Fs.e^m < 4. Given the uncertainties, their IRAC 
colors are equally consistent with elliptical galaxies at 
z > 2 or moderately obscured (Ay « 2.5) z > 1 star- 
bursts. However, all possess radio-loud AGN by virtue 
of their negative (/-values and so cannot be pure star- 
bursts. 

Figure 8 highlights the fact that the IRAC detected 
OIRSs as a group have much flatter mid-IR SEDs than 
the OIMSs. The latter possess mean F24/im/-p3.6^m ra- 
tios of 170, 65, 52, and 44 for the obsSy2, obsSyl, sb, 
and comb classes respectively. By contrast, the thirteen 
IRAC detected OIRSs without z p hot (i.e., the filled circles 
in Figure 8) possess a mean F 2 4 Mm / F 3 , 6lim < 11. Given 
their 3.6 /im flux densities, all thirteen should have been 
detected at 24 /tm had their SEDs been steep power-laws 
like the obsSyl and obsSy2 OIMSs. 7 Even if their SEDs 
resembled the sb or comb OIMSs, most would have been 
detected at 24 /im. This is consistent with emission from 
stellar photospheres and AGN heated dust being at least 
comparable in most of the IRAC detected OIRSs. 

3.3. Luminosities and Stellar Masses 

Rest-frame K-band luminosities, defined as Lk = 
Ai/L v (2.2/U.m), with Av corresponding to the Johnson 
K-band filter FWHM (2.1 x 10 13 Hz; Tokunaga et al. 
2002) were calculated for the OIMSs and the six OIRSs 
with photometric redshifts. These are listed in Tables 3 
and 4 in units of the sun's K-band luminosity (L Qi k = 8.7 
x 10 24 W; Bessel 1979). With this definition, Arp 220 
and IRAS F00183-7111 would have Johnson K-band lu- 
minosities of 1.3 x 10 11 L Qi k and 9.1 x 10 11 L ,k respec- 
tively. Monochromatic rest-frame luminosities at 5 /im, 
defined as L 5Mm = j/L I/ (5/im), were also calculated for 
these sources. This quantity provides a measure of the 
emission from hot (400 — 1500 K) dust free of any con- 
tribution from evolved stars and which avoids the (rest- 
frame) 9.7/irn silicate feature. These are also given in 
Tables 3 and 4. Uncertainties in redshift constitute the 
dominant source of error in the derived luminosities. 

The OIMSs arc highly luminous in the rest-frame K- 
band, with a median Lk of 1.6 x 10 12 L 0i k and more 
than an order of magnitude spread. The obsSy2 OIMSs 
are the most luminous, with Lk ~ 10 13 L 0i k on average. 
These are followed by the obsSyl and comb OIMS, both 
with average Lk ~ 10 12 L Qi k- The obsSy2 and comb 
OIMSs are both on average an order of magnitude more 
luminous than their local templates, IRAS F00183-7111 

7 For example, OIRS #232, with the smallest F 3 . 6Mm = 8.1 ± 
1.7 ^iJy, would have yielded 24 (im flux densities between 0.5 and 
1.4 mjy, easily detectable by the MIPS survey. 
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and Arp 220. Their rest-frame monochromatic 5 fim lu- 
minosities are also extremely large, with median L5 Mm = 
1.0 x 10 39 W and 1.2 x 10 38 W for the obsSy2 and obsSyl 
OIMSs respectively. 

The six IRAC detected OIRSs with photometric red- 
shifts show a very wide range in Lk, led by #176 with 
2.8 x 10 13 Lq,k- This is roughly fifty times more lumi- 
nous than the five other OIRSs and nearly twice as large 
as the most luminous OIMS. The three OIRSs that were 
best fit with a comb template have a mean Lk of 5 x 10 11 
L Q! k- This is within a factor of two of the average for the 
three comb OIMSs (OIMS #7, 12 and 14). Accurate 8- 
1000 (im luminosities (Lir) for the OIMSs and OIRSs are 
not possible since none were detected at 70 /im or 160 
/im. However, the fact that median Lk for the OIRSs 
and comb/sb OIMSs are similar suggests that they have 
comparable total luminosities, and that at least some of 
the IRAC detected OIRSs are hyper-luminous. 

The 1.6 /im bump apparent in the averaged OIRS SED 
in Figure 7 points to a significant evolved stellar compo- 
nent. For these objects we estimated stellar masses by 
multiplying their Lk by near-infrared mass to light ra- 
tios derived from burst models in Bell & dc Jong (2001). 
From their Figure 1, we adopted M*/Lk = 0.3, which 
corresponds to their bluest models (B — R = 0.6). The 
mass to light ratio increases fairly slowly with color, as 
can be seen from the fact that M*/Lk = 0.5 at B — R 
= 1.0. We find a wide range in derived stellar masses, 
again led by OIRS #176 with M* w 10 13 M . The other 
OIRSs have masses in the range 0.8 — 3.4 x 10 11 M , with 
a mean of«2xl0 11 M Q . Very massive galaxies are in- 
dicated even if warm dust makes a 50% contribution to 
the rest-frame 2.2 /zm emission. 

Despite the lack of redshifts for the remaining IRAC 
detected OIRSs, we can set useful limits to their stellar 
masses if we assume they lie at z > 1 and that the (pu- 
tative) AGN makes an insignificant contribution to the 
IRAC 3.6 //m and 4.5 fim bands. The first assumption 
is reasonable given their faint optical limits. The second 
follows from the data in Table 1, where only OIRSs #49, 
79, 232, & 410 have mid-IR to radio SEDs that might al- 
low significant synchrotron emission in the mid-IR. For 
z = 1, where the IRAC 4.5 band measures rest-frame 
2.2 fim emission, the average i<4.5^ m (18.0 /zJy) implies a 
K-band luminosity of 1.5 x 10 11 Lq.k- For an early-type 
M*/Lk ratio, we derive a typical stellar mass of « 5 x 10 10 
M Q for this sub-population. So long as there is minimal 
contribution from an AGN to the 4.5 /im emission, the 
IRAC detected OIRSs without redshifts also appear to 
be massive galaxies. 

For the majority of OIMSs, the rest-frame near- 
infrared is dominated by non-stellar emission, and the 
upper-limits to their stellar masses are not very infor- 
mative. However, the four OIMSs with comb or sb IRS 
spectra (#2, 7, 12, and 14) show evidence of a weak 1.6 
/zm bump in Figures 1 and 4. Using the same starburst 
M»/Lk of 0.3, their rest-frame Lk imply maximum M* 
between 2 — 4 x 10 11 M Q . These objects also appear to 
be very massive even if the stellar component makes only 
a ~ 20% contribution to the rest-frame near-infrared. 

Rest-frame 20cm luminosity densities were estimated 
for the OIMSs and the six IRAC detected OIRSs with 



Zp ho t using 

L 20 cm = 1-2 x 10 20 (1 + z)"" 1 F„(20cm) W Hz" 1 , 

(3) 

where Dl is the luminosity distance in Mpc, a is the ra- 
dio spectral index (F v cx z^~ Q ), taken to be 0.65, and 
F u (20cm) is the observed 20 cm flux density (or 5<r 
upper-limit) in Jy. The latter are taken from Higdon05 
and de Vries et al. (2002). Seventy per cent of the OIMSs 
were not detected in the radio, implying L 2 o cm upper- 
limits of - 5 x 10 24 W Hz" 1 (5a). Of the remainder, 
OIMS #13 in Table 4 is noteworthy for possessing the 
largest radio luminosity density among both OIRS and 
OIMS, with L 20 cm = l.Ox 10 26 W Hz" 1 . It represents the 
sole OIMS that is squarely in the "powerful" radio galaxy 
class (i.e., L 20 cm > 10 25 W Hz" 1 ; McCarthy 1994). A 
wide range of values is found for the OIRSs with photo- 
metric redshifts, with L 20c m between 0.2-2.9 x 10 25 W 
Hz" 1 and a median of 4 x 10 24 W Hz" 1 . 

The six IRAC detected OIRSs show signs of large 
evolved stellar populations by virtue of the 1.6 ^m bump 
apparent in Figure 7. If we assume that their radio emis- 
sion arises primarily from supernova remnants, we can 
estimate their current star formation rates with their 
L 20c m- Excluding the radio-loud OIRS #97 (q = -0.24), 
we derive a median SFR of 2200 M© yr" 1 after Bell 
(2003). For the single sfc-like OIMS the L 20cm upper- 
limit implies SFR < 2000 M Q yr _1 . Assuming a mini- 
mal contribution to L 20cm from an active nucleus, we find 
these objects to possess SFRs comparable to the highest 
values derived for local ultra-luminous infrared galaxies 
(ULIRGs) using both X-ray and far-infrared luminosities 
(e.g., Persic et al. 2004). 

3.4. The Near Environments of OIMSs and OIRSs 

Most of the IRAC detected OIRSs appear to possess 
nearby companion candidates, whereas few of the OIMSs 
do. Figure 3 shows that the OIMSs are typically unre- 
solved by IRAC at 3.6 fim, with only OIMS #3 and 
6 showing hints of structure. By contrast, twelve of the 
nineteen IRAC detected OIRSs (63±18%) in Figure 1 are 
clearly not simple point-sources. Examples include OIRS 
#79, 97, and 363. While gravitational lensing by fore- 
ground clusters can not be ruled out, the most straight- 
forward interpretation is that IRAC is detecting several 
closely spaced galaxies. This is consistent with the ma- 
jority of the IRAC detected OIRSs having companions 
within w 50 kpc. 

The clustering properties of these systems are beyond 
the scope of this paper. It is worth noting, however, that 
neither OIMSs nor OIRSs appear to inhabit particularly 
crowded environments like the cores of rich clusters. For 
example, using the sixteen IRAC subimages of OIMSs 
in Figure 3 we find an average of 1.31 ± 0.29 sources 
per 20" wide field (not including the OIMSs). The aver- 
age number of sources in the fields centered on OIRSs is 
1.74 ± 0.30, again not including the OIRSs. As a base- 
line figure, we would have expected 1.34 ± 0.01 sources 
in a random 20" sized region given the total number of 
sources (« 370,000 at 3.6 ^m) detected within the 8.5 
deg 2 IRAC Shallow Survey (Eisenhardt et al. 2004). 8 

8 We used the Student t-Test to evaluate the likelihood that 
the IRAC detected OIRSs have a significantly higher number of 
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This is not to suggest that OIMSs and the IRAC de- 
tected OIRSs do not inhabit clusters. Recent work indi- 
cates that radio galaxies and quasars show a wide variety 
in local density up to at least z « 1.5, from the cores of 
rich clusters to field- like environments (e.g., Barr et al. 
2004; Best et al. 2003). Further, the radio galaxies and 
radio-loud quasars that are associated with rich centrally 
concentrated clusters are as often as not situated outside 
the actual core. Given these trends and the rather small 
spatial dimensions of the IRAC subimages (area « 0.04 
Mpc 2 ) it is perhaps not surprising that we do not see an 
excessive number of companion galaxies surrounding the 
OIMSs or OIRSs. An investigation of the environments 
of OIMSs and OIRSs on scales of several Mpc 2 will be 
presented in a future paper. 

4. DISCUSSION 
4.1. OIMSs as Buried QSOs 

Houck05 concluded that both the obsSyl and obsSy2 
OIMSs are heavily obscured and hyper- luminous, with 
Lir estimated to be 10 13 Lq. Their IRS spectral 
shapes, and in particular, the absence of PAHs, im- 
plied that their primary energy source was an active nu- 
cleus rather than star formation. We have shown that 
the obsSyl and obsSy2 OIMSs also possess steep power- 
law SEDs over the rest-frame 1 — 10 /im region and ex- 
tremely high rest-frame monochromatic 5 /im luminosi- 
ties (L 5([lm = 10 38 -10 39 2 W). The latter originates in the 
hottest dust component and necessarily implies a power- 
ful active nucleus. To place these luminosities in context, 
we calculated rest- frame isL u (5(j,m) for the sample of 42 
3CRR radio galaxies and quasars observed with Spitzer 
by Ogle et al. (2006). 9 For the sources classified as 
mid-IR luminous, which consist of quasars and narrow- 
line radio galaxies (NLRGs), L 5Aim = 10 374 - 10 39 3 W. 
The quasar-like L5 Mm of the obsSyl and obsSy2 OIMSs, 
together with their steep power-law SEDs and IRS spec- 
tra provide compelling evidence that they are driven by 
powerful AGN accreting near the Eddington limit. These 
sources may be considered to be obscured or "buried" 
quasi-stellar objects (QSOs). 

Final proof that the obsSyl and obsSy2 OIMSs host 
heavily obscured QSOs would be the detection of lumi- 
nous hard (i.e., > 2 kcV) X-ray emission. These pho- 
tons originate very near the super-massive black hole and 
are a reliable indicator of powerful AGN activity. They 
are also fairly insensitive to obscuration, at least up to 
Nh = 1.5 x 10 24 cm -2 where the optical depth for Comp- 
ton scattering reaches unity. Of the OIMSs discussed in 
this paper, only #13 (not coincidently, the most radio- 
loud OIMS) was detected above the 99% confidence level 
in the "hard" (F 2 _ 7 kcV = (1.6 ± 0.8) x 10~ 14 erg s" 1 
cm" 2 ) and "full" (F .5-7 kcV = (2.3±1.4) x 10~ 14 erg s" 1 
cm~ 2 ) bands in the 2 ks Chandra survey of the Bootes 

companion candidates within the 20" wide fields relative to the 
OIMSs. For the sources visible in the two sets of IRAC sub-images, 
we calculated t = 1.43. For thirty-two degrees of freedom and a 95% 
confidence level, t would have had to exceed 2.04 for a significant 
difference between the OIMSs and OIRSs near environments to be 
real. Wc thus find no measurable difference in the apparent density 
on this scale. 

9 Rest-frame ^L„(5/im) were extrapolated for the quasars and 
radio galaxies in Table 1 of Ogle ct al. using the listed rest-frame 
vLv(15fim) and the 7 to 15 (im spectral indices. 



field by Renter et al. (2005). This corresponds to a 
hard X-ray luminosity L2-7 kcV = (5±2) x 10 37 W, with 
no correction for intrinsic extinction. 10 For the other 
eight obsSyl and obsSy2 OIMSs we obtain upper-limits 
of L2-7 kcV i$ 10 37 - 5 W, again with no extinction cor- 
rection. Powerful AGN in the local universe typically 
possess L2-7 kov « 10 35 - 10 37 W. So while only OIMS 
#13 can be conclusively shown to be a "buried" QSO 
based on its hard X-ray emission, the other obsSyl and 
obsSy2 OIMSs are still consistent with this interpreta- 
tion. Deeper observations (or a stacking analysis) to de- 
tect hard X-rays with Chandra should be carried out, as 
well as searches for high-ionization infrared lines such as 
[S IV] A10.51 fxm and [Ne V] A24.31 fxm. These were not 
evident in the IRS spectra shown in Houck05, though 
this can be attributed to the low spectral resolution cou- 
pled with the relative faintness (F 2 4 Mm w 1 mJy) of the 
sources. 

4.2. The OIRSs Detected by IRAC 

The IRAC detected OIRSs fall into two groups depend- 
ing on whether or not they were also detected with MIPS 
at 24 /im. Those that are detected at 24 (im appear to 
be massive and dusty z « 2 galaxies that are primar- 
ily driven by star formation. Those not detected at 24 
/im (by far the majority) possess radio-loud AGN as ev- 
idenced by their negative q upper-limits. They also ap- 
pear to be massive, with M* w 10 11 M Q if they are at 
z > 1 as their faint optical limits imply. Their individual 
-p24/xm/ F3.6fj.rn limits are consistent with mid-IR SEDs in 
which stellar emission is at least comparable to that from 
dust close to the central engine. 

Figure 9 shows the ratio i 7 24 ( um/-p3.6^m as a function of 
rcdshift for five representative objects known to contain 
powerful AGN: the composite obsSy2 OIMS from Figure 
7, the NLRG 3C459, the quasar 3C273, the Low Exci- 
tation Radio Galaxy (LERG) 3C293, and the "red and 
dead" radio galaxy LBDS 53W091 (Stern et al. 2006). 
Except for the OIMS and LBDS 53W091 the SEDs were 
constructed using data from NED. This figure also il- 
lustrates the extremely steep A rost = 1 — 10 /im SED 
of the obsSy2 OIMSs relative to a fairly typical NLRG 
and quasar. The thick horizontal line in the figure repre- 
sents the mean F24 fim / F 3 , &flm upper- limit for the IRAC 
detected OIRSs (i*24/xm/-F3.6Aim < 3.9). This was deter- 
mined by dividing the stack-averaged 24 /im upper-limit 
(67 /iJy, 3c) by the mean F^q^ for the twelve OIRSs not 
detected at 24 /im. The upper/lower bounds to these lim- 
its, defined by dividing the stack-averaged 24 /im limit by 
the minimum/maximum F3. Q^ m values are also shown as 
two thin horizontal lines. Taken as a group, these OIRSs 
have spectral shapes that are inconsistent with obsSy2 
OIMSs, typical NLRGs, and quasars at essentially all 
rcdshifts. They are, however, consistent with LERGs 
and the "red and dead" radio source LBDS 53W091 for 
z>l. 

Of the two, LERGs arc currently the best studied at 
infrared wavelengths. LERGs are identified by their nu- 
clear spectra which show extremely weak high-excitation 
optical lines like [O III] A5007 A (Laing et al. 1994), and 

10 X-ray luminosities were calculated from L2— 7 keV = 
47rD^ -F2-7 kev(l + 2 ) r ~ 2 , where Dl is the luminosity distance 
and r is the power-law index, taken to be 1.8. 



8 



Higdon et al. 



are nearly always hosted by massive ellipticals. Most FRI 
radio galaxies are LERGs, though a substantial number 
of FRIIs are as well. LERGs display a wide range of ra- 
dio luminosities, with L 20c m ~ 10 23 W Hz" 1 to > 10 26 
W Hz -1 , showing that LERGs can be "powerful" ra- 
dio galaxies. Their flat rest-frame near to mid-IR SEDs 
have been attributed to the absence of massive obscur- 
ing tori like those inferred in NLRGs and BLRGs. More 
direct evidence comes from X-ray spectra of LERG nu- 
clei, which (with few exceptions) lack an absorbed soft 
component (Evans et al. 2006), and radio VLBI observa- 
tions showing polarized emission on parsec scales (Kharb 
et al. 2005). 11 Sicbcnmorgen et al. (2004) used simple 
dust models to fit the optical/infrared SEDs of a number 
of LERGs in the 3CRR (Laing et al. 1983). They con- 
cluded that compared with NLRGs and BLRGs, LERGs 
possess lower AGN luminosities, more extensive dust dis- 
tributions (i.e., much less dust at small radii), and lower 
overall obscuration (Ay ~ 2). 

The small F 2 4 Mm /F 3 . 6Mm limits for the OIRSs not de- 
tected at 24 /zm are consistent with a population of high- 
z LERGs, i.e., massive elliptical galaxies hosting radio- 
loud AGN with relatively small amounts of dust close to 
the central engine. However, the data are also consistent 
with a population of radio-loud z > 1 ellipticals with ac- 
tive nuclei weak enough to be "drowned out" by stellar 
emission, even if only temporarily, due to a drop in the 
mass accretion rate. 

Whether the OIRSs discussed in this section are 
LERGs or not, their radio luminosities will be substan- 
tially less than 10 25 W Hz" 1 if their typical redshifts are 
w 1. Only if z > 2 will most of these objects be classi- 
fied as "powerful" radio galaxies. Future observations to 
better constrain their rest-frame infrared SEDs and de- 
termine accurate redshifts will be required to help resolve 
this issue. 

It may at first seem puzzling that so many of the IRAC 
detected OIRSs resemble LERGs and "red and dead" ra- 
dio galaxies. One may ask where are the optically "invis- 
ible" mid-IR luminous NLRGs? Part of the answer must 
lie in the way the OIRS sample was originally defined, 
namely, by requiring a compact or unresolved radio mor- 
phology in the A-array 20 cm images. The original in- 
tent of the VLA survey was to detect optically obscured 
z > 1 starburst galaxies through their radio emission. 
Such sources were expected to be point-like in the 20 cm 
A-array images. As a result, extended radio structures 
(i.e., jets and lobes) were excluded from consideration 
from the outset since they (1) would identify AGN dom- 
inated sources, and (2) would have also proven difficult 
to confidently associate with an optical counterpart given 
the high source density in the NDWFS B\y, R, and I- 
band images. As a consequence, it is possible that we 
preferentially selected beamed radio sources, i.e., sources 
resembling BL Lac objects. It has been suggested that 

11 This has led to suggestions that the accretion process in 
LERGs is fundamentally different from NLRGs and BLRGs (e.g., 
Hardcastlc et al. 2006). Accordingly, super-massive black-holes 
in NLRGs and BLRGs accrete material from a dusty molecular 
torus while those in LERGs accrete hot gas directly from the in- 
tcrgalactic medium through Bondi accretion (Allen et al. 2006). 
For our purposes, this is less important in understanding their in- 
frared SEDs than the apparent lack of a massive obscuring torus 
extending close to the central engine. 



FRI radio galaxies are the parent population for BL Lac 
objects. The fact that LERGs comprise most of the lo- 
cal FRI allows a simple explanation, namely, that most 
of the IRAC selected OIRSs are in fact a population of 
high-z BL Lac objects. 

4.3. The OIRSs Not Detected by IRAC 

The fact that nearly half (16/35) of the OIRSs appear 
to lack IRAC counterparts (see Figure 2) is a significant 
result, especially in light of the small number of spurious 
radio sources expected in Higdon05. 12 By contrast, all 
of the OIMSs were detected by IRAC, which suggests 
that the two groups have very different SEDs and/or 
redshifts. We can constrain the average Lk of the non- 
IRAC detected OIRSs (and indirectly, their likely red- 
shift range) with the stack-averaged 3er upper-limits at 
3.6, 4.5, 5.8, and 8.0 fim listed in §2. 

If the OIRSs not detected by IRAC are primarily a 
z < 1 population, they must be sub-L* galaxies. For 
an object at z — 0.64, where rest-frame 2.2 ^m shifts 
to the center of the IRAC 3.6 ^m band, the 1.5 /iJy 
upper-limit would correspond to Lk < 1-2 x 10 10 L 0j k- 
This is substantially smaller than the K-band luminosi- 
ties of L* galaxies determined by Kochanek et al. (2001), 
who found L^ = 0.7-1.2 x 10 11 L QjK for early to late 
types. 13 This in turn limits the stellar mass of the non- 
IRAC detected OIRS hosts. For an early-type galaxy 
M*/L K « 0.7 (Bell & de Jong 2001), and the resulting 
M* is less than 8 x 10 9 M . Still smaller masses result 
for later host types since M*/Lk will only decrease for 
these galaxies. Galaxies of this mass are unlikely to host 
radio- loud AGN as powerful as « 10 24 W Hz" 1 , which 
is the inferred average L2o C m for the non-IRAC detected 
OIRSs over 0.5 < z < 1.0. 

On the other hand, if the host galaxies of these OIRSs 
are z « 1 — 2 ellipticals as suggested by Gruppioni et 
al. (2001), we would have expected to detect most with 
IRAC. A spinning super-massive black hole is generally 
held to be a requirement for a radio-loud active nucleus 
(Rees 1984). Since the black hole mass correlates with 
the mass of the host galaxy's stellar bulge (e.g., Magor- 
rian et al. 1998), radio-loud AGN hosts are known to 
be quite massive, and considerably more luminous than 
L*. Using the above stack-averaged limits, we find the 
rcdshift where a massive elliptical host galaxy with L K 
= 3L K would escape detection to be z > 2, with a corre- 
sponding average L 2 o C m of « 10 25 W Hz" 1 . The OIRSs 
not detected by IRAC appear to represent a population 
of powerful radio galaxies at z > 2. 

12 There are 392 radio sources in the final catalog. The number 
of false radio sources was estimated by repeating the SExtractor 
runs on the final VLA images multiplied by -1. The number of > 5<r 
negative "sources" provides a robust estimate of the number of false 
positives. Only three such sources were found in the radio images, 
a number that is also consistent with random fluctuations given 
the number of independent synthesized beams in the radio images. 
Minimal contamination of the source catalog is indicated. We also 
considered the possibility that some of the OIRSs might be hot- 
spots in large radio doubles "resolved out" by the interferometer 
array. However, given their wide separation (few OIRSs arc within 
2-3 minutes of each other) these would have to be gigantic radio 
sources, and therefore unlikely. 

3 Kochanek et al. derived using a complete sample of 4,353 
galaxies in the 2MASS survey. For H = 71 km s _1 Mpc — 1 , their 
results translate into MJ^ = -23.7 (late-type), -24.3 (ear ly- type) , 
and -24.1 (all types). 
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This may not be the entire story, of course. The OIRSs 
lacking IRAC counterparts show a substantially wider 
dispersion in i^ocm compared with the OIRSs detected 
by IRAC. Of the four OIRSs with F 20c m > 1 mJy, three 
lack IRAC counterparts. At the same time, eight of the 
thirteen OIRSs fainter than 0.25 mJy at 20cm lack IRAC 
counterparts. It may well be that the non-IRAC detected 
OIRSs consist of two distinct populations: one charac- 
terized by L 2 ocm < 10 25 W Hz -1 and z w 1 — 2 (corre- 
sponding to the i*20cm < 0.25 mJy objects) and a much 
more powerful and distant sub-group of radio galaxies 
(i.e., those with i^ocm > 1 rnJy). Further progress in 
understanding the nature of these objects will require ei- 
ther photometric redshifts derived with deeper Spitzer 
IRAC and MIPS observations, or through ground-based 
spectroscopy using the next generation of large optical 
telescopes. 

4.4. Comparison with Sub- Millimeter Galaxies 

It is worth comparing properties of the OIMSs and 
OIRSs with other luminous objects at z « 1 — 4, in par- 
ticular, the galaxies discovered at sub-millimeter wave- 
lengths using instruments like SCUBA (Smail et al. 
1997). Sub- millimeter galaxies (SMGs) are extremely lu- 
minous (Lir w 10 13 L ) and dusty systems, with a red- 
shift distribution that peaks near z — 2.5 (Chapman et 
al. 2003). Most appear to be powered primarily by star 
formation. Indeed, a significant fraction of star formation 
in the early universe may take place in them. SMGs also 
tend to be extremely faint or "invisible" at optical wave- 
lengths (e.g., Hughes et al. 1998). To ensure meaningful 
comparisons with the OIMSs and OIRSs, we selected a 
sample of ten SMGs detected at 850 /xm with SCUBA in 
the Lockman Hole East region that were observed with 
both IRAC and MIPS by Egami et al. (2004). Four of 
the SMGs were also detected at 1.2 mm with MAMBO 
(Eales et al. 2003) by Ivison et al. (2004), and nine 
were detected at 20 cm using the VLA's A-configuration 
(Ivison et al. 2002). The rcdshift range for the SMG 
sample (z « 1 — 3) is similar to that of the OIMSs and 
the IRAC detected OIRSs with z pno t- Rest- frame radio 
and infrared luminosities as defined in §3.3 were calcu- 
lated for the SMGs and are listed in Table 5. Likewise, 
averaged rest-frame 1 — 10/xm SEDs were calculated as 
in §3.2, and are shown in Figure 7 for the "cold" SMGs. 

Two of the ten SMGs possess power-law SEDs in the 
IRAC and MIPS data, suggesting a significant AGN con- 
tribution ("warm" sources in Egami et al.'s terminol- 
ogy). Their power-law exponents (a — —1.9) are sim- 
ilar to those of the obsSyl OIMSs. The largest K-band 
luminosity belongs to the "warm" SMG LH850.8b (L K 
= 5.3 x 10 12 L q ,k), which is within a factor of » 2 of 
the median Lk for the hyper-luminous obsSy2 OIMSs 
(f=s 10 13 L 0iK ). The remaining eight SMGs have SEDs 
consistent with obscured star formation providing the 
bulk of their luminosity (i.e., "cold" sources in Egami et 
al.'s terminology). 

There is considerable overlap between the emission 
properties of the "cold" SMGs and the sb/comb OIMSs 
and IRAC detected OIRSs with z p hot- For all three 
groups we derive comparable luminosities at rest-frame 
K-band (average Lk = 8 x 10 11 L 0! k for the "cold" 
SMGs) and 5 ^m (average i/L„(5 /jm) « 10 38 W), both 
over an order of magnitude smaller than the obsSyl and 



obsSy2 OIMSs. The similarity between the "cold" SMGs 
and the sb/comb and IRAC detected OIRSs with z p hot 
extends to their averaged rest-frame 1 — IO^to SEDs as 
well, as Figure 7 shows. In particular, both the "cold" 
SMGs and IRAC detected OIRSs with z pno t show a 
prominent 1.6 fim bump implying substantial evolved 
stellar components. This is followed by a gradual rise to 
A rcst ~10 fim as expected for dusty starbursts. Adopting 
a starburst M*/Lk of 0.3, we estimate a median stellar 
mass of M* = 2 x 10 n M Q for the "cold" SMGs, which 
is comparable with the starburst dominated OIMSs and 
OIRSs. Similarly assuming that the radio continuum 
from the "cold" SMGs originates in supernova remnants, 
star formation rates between 500 and 1700 M Q yr _1 (me- 
dian SFR = 1100 M Q yr _1 ) are implied by their L 2 o C m- 
The similar luminosities, implied star formation rates 
and stellar masses, faint optical magnitudes, and SED 
shapes is evidence that at least a fraction of the OIMSs 
and IRAC detected OIRSs are members of the same par- 
ent population as the "cold" (i.e., starburst dominated) 
SMGs. 

The top panel of Figure 8 includes the ten "warm" 
and "cold" SMGs in the IRAC two-color diagram, along 
with the OIMSs. As expected, the two "warm" SMGs 
are positioned in the region occupied by the AGN pow- 
ered obsSyl OIMSs. Neither have mid-IR SEDs that 
rise as fast as the obsSy2 OIMSs, though with only two 
"warm" sub-millimeter sources it is impossible to draw 
any definite conclusions. The eight "cold" sources are 
situated close to the dusty M 82 track at z ~ 1 — 3, 
both consistent with their redshifts and with obscured 
star formation providing the bulk of their luminosities. 
There is a clear separation between the "cold" SMGs on 
one hand and the AGN dominated obsSyl and obsSy2 
OIMSs on the other. 

Figure 10 shows the relative positions of the OIMSs, 
OIRSs, and the ten Lockman Hole East SMGs in a ra- 
dio/infrared two-color diagram. Systematic trends can 
be seen: (1) obsSyl and obsSy2 OIMSs extend to the 
right of center (F%. 0)lra / F 3 . &flm > 3.5) due to their steep 
power-law SEDs; (2) comb and sb OIMSs are situated 
primarily to the left of this, with Fs.Ofim/ F3.6fj.u1 ratios 
between 1.6 and 4; (3) "cold" SMGs are left of center 
(Lg.o^m/Ls.e^m < 3); (4) "warm" SMGs are to the right 
of the "cold" ones (F 8 . 0Aim /F 3 . 6Aim > 4); (5) the IRAC 
detected OIRSs primarily occupy the lower-left portion 
of the diagram. A few of the latter overlap the region 
occupied by the "cold" SMGs. However, the majority of 
IRAC detected OIRSs are radio-loud (13/19 have q < 0) 
and are well separated from five of the six OIRSs with 
•Zphot- It is also apparent that with two exceptions, the 
OIMSs and the "warm" and "cold" SMGs obey the radio- 
infrared correlation within the uncertainties. 1 OIMSs, 
OIRSs, and SMGs thus occupy distinct regions in a ra- 
dio/infrared two-color diagram, though there is also con- 
siderable overlap, particularly between the OIMSs and 
SMGs. 

This overlap is further illustrated in Figure 11, where 
Lk is plotted against i*8.0jum/-F3.6^m- A clear progression 
is seen, with luminous and flat-spectrum "cold" SMGs 

14 OIMS #13 is the most radio-loud OIMSs in Houck05, while 
the "cold" SMG [LE 850] 8a may have strong PAH emission within 
the 24 (im passband, giving rise to its large Rvalue. 
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(i.e., L K « 10 11 - 10 12 L QiK and F 8 . 0iUm /F 3 . 6iUm ~ 1 - 2) 
giving way to obsSyl OIMSs and "warm" SMGs (Lk ~ 
few x 10 12 L , K and F 8 . ^ m / 'F^.q^ ~ 5-10), and fi- 
nally to the steep power-law obsSy2 OIMSs (L K > 10 13 
Lq,k and Fs.o^m/ F 3 . 6fim > 10). This suggests that the 
OIMSs, a subset of the IRAC detected OIRSs, and the 
SMGs are all taken from the same population of dusty 
hyper-luminous high-z objects. The chief difference be- 
tween them is the dominant emission mechanism: obsSyl 
and obsSy2 OIMSs in particular can be thought of as 
analogs of SMGs sources whose luminosities in the rest- 
frame near and mid-IR are dominated by powerful and 
heavily obscured active nuclei accreting near the Edding- 
ton limit. 

Figure 10 also illustrates a dichotomy between the 
IRAC detected OIRSs with and without rcdshifts (i.e., 
those with and without 24 /im detections). With one 
exception, all of the former are characterized by q > 0, 
which is consistent with systems dominated by star for- 
mation, a radio-quiet AGN, or a combination of the two. 
By contrast, the remaining IRAC detected OIRSs all 
have negative ^values indicating an increasingly dom- 
inant radio-loud active nucleus. The parent OIRS popu- 
lation is similarly dominated by radio-loud active nuclei, 
for excluding the six OIRSs with %,hot> 72% of the sources 
in Table 1 of Higdon05 have negative q. This suggests 
that the 24 /xm detected OIRSs in Figure 5 are not rep- 
resentative of the larger OIRSs population, but rather 
form a starburst dominated sub-group. 

4.5. Do OIMSs and SMGs Form an Evolutionary 
Sequence ? 

Results from the previous sections suggest that the ob- 
sSyl/ obsSy2 OIMSs and "cold" SMGs arc members of 
the same parent population of heavily obscured hyper- 
luminous galaxies, with the fundamental difference be- 
ing their dominant power source, i.e., either accretion 
near the Eddington limit or star formation. The ob- 
sSyl and obsSy2 OIMSs in particular can be thought 
of as heavily obscured QSOs. While it is true that very 
few SMGs contain "buried" QSOs (fa 5%, Almaini ct 
al. 2003), recent ultra-deep X-ray observations argue 
persuasively that many SMGs contain massive central 
black holes and luminous, though not yet dominant, ac- 
tive nuclei. For example, 75% of the SMGs with radio 
counterparts and known rcdshifts in the Chandra Deep 
Field North were detected in hard X-rays by Alexander 
et al. (2005). Their X-ray luminosities (Lo.5-8 kcV ~ 
10 36 — 10 37 W), levels of obscuration (eighty percent 
show Nh ~ 10 20 — 10 24 cm -2 ), and spectral power- 
law exponents (r fa 1.8) are similar to nearby powerful 
AGN. However, their X-ray to infrared luminosity ratios 
(Lx/Lir fa 0.4%) are typically an order of magnitude 
smaller than QSOs, consistent with intense star forma- 
tion (SFR~ 1000 M yr -1 ) providing the bulk of their 
luminosity. This may be simply saying that that the X- 
ray detected SMGs are currently accreting at substan- 
tially sub-Eddington rates, which would imply that the 
only real difference between the obsSyl/ obsSy 2 OIMSs 
and "cold" SMGs is the rate at which material is trans- 
fered to the central engine. On the other hand, near-IR 
spectra of SMGs typically show that when broad rest- 
frame optical lines are present they are typically only 
fa 1000 - 3000 km s" 1 wide (cf. Swinbank et al. 2004), 



which is consistent with only modest black hole masses, 
i.e., < 10 8 M Q . It is thus possible (though admittedly 
speculative) that what distinguishes the AGN powered 
OIMSs and "cold" SMGs is primarily the central black 
hole mass rather than simply the accretion rate. 

Given the well established correlation between the 
masses of the spheroid and central black hole in local 
galaxies (e.g., Kormcndy & Richstone 1995; Magorrian 
et al. 1998) and the evidence linking the SMG population 
with current epoch massive ellipticals (e.g., Smail, Ivison, 
& Blain 1997; Barger et al. 1998; Barger et al. 1999), it 
is intriguing to consider the possibility that the obsSyl 
and obsSy2 OIMSs represent SMGs that have made the 
transition from a starburst dominated "cold" phase to an 
accretion dominated "buried" QSO phase. A plausible 
mechanism for this transition would be the growth of a 
sufficiently massive black hole (in parallel with the stellar 
bulge) during the high-SFR "cold" phase, with the subse- 
quent transfer of material to the central regions through 
tidal interactions or mergers. The energy released dur- 
ing this relatively brief phase would act to expell the ob- 
scuring gas and dust from the central regions, effectively 
halting the growth of the bulge and super-massive black 
hole and "unveiling" a QSO in the manner envisoned by 
Sanders et al. (1988). 

SMGs and OIMSs appear to have have very differ- 
ent mean space densities. If so, this would have im- 
portant implications for the envisoned SMG-QSO tran- 
sition. For example, the number of SMGs per square 
degree with F 8 5o Mm > 5 mJy (which we take to represent 
the hyper- luminous sources) is w400 using the cumula- 
tive 850 fxm source count formulation in Sanders (2000). 
Assuming they all lie between z fa 1 — 3, their average 
space density becomes psmg = 2 x 10~ 5 Mpc~ 3 . There 
are seventy sources in the NDWFS Bootes field satisfy- 
ing F 2 4^ m > 0.75 mJy and R > 24.5, i.e., OIMSs by 
our definition, of which fifty-eight have been observed 
by Spitzer. Forty-five of these OIMSs can be confidently 
classified as obsSyl or obssSy2 on the basis of their IRS 
spectra and/or steep power-law infrared SEDs (Higdon 
et al. in preparation). Over half of these sources have 
z S pcc = 0.8 — 3.1 on the basis of IRS observations, and 
the remainder have IRAC/MIPS SEDs consistent with 
this redshift range. 15 There are thus « 50 obsSyl 
and obsSy2 OIMSs in the 8.5 deg 2 NDWFS. Assum- 
ing that the OIMSs with no redshift measurements are 
within z fa 1 — 3, we estimate a mean space density of 
Poims ~ 3 x 10~ 7 Mpc -3 for this population, or nearly 
two orders of magnitude smaller than psmg- This could 
be simply explained if the duration of the "buried" QSO 
phase is short compared with the SMG's "cold" phase, 
i.e., the time spent growing a sufficiently massive black 
hole and bulge. 

One caveat to the above analysis is the sensitivity of 
Psmg to the adopted 850 /im flux density cut-off. For ex- 
ample, if only the Fsso/xm > 10 mJy SMGs pass through 
a obsSyl or obsSy2 phase, then psmg ~ 5 x 10~ 6 Mpc -3 , 
and Psmg/poims ~ 17- The sub-millimeter properties of 
the OIMS population are largely unknown, and it is not 

15 It was noted in §3.2 that the obsSy2 OIMSs show a pro- 
nounced break in their power-law SEDs near A res t ~ 5 — 8fim. For 
the obsSy2-\ikc OIMSs in the larger sample with featureless IRS 
spectra, the observed wavelength of this SED break was found to 
be consistent with z ~ 1 — 3. 
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certain whether OIMSs with i^^m > 0.75 mJy map onto 
SMGs with F 8 5o^ m > 10 mJy, or for that matter, SMGs 
with F85o M m > 2 mJy. 16 Despite these uncertainties, it 
appears likely that Psmg/poims ~ 20 — 70. This suggests 
that if indeed the obsSyl and obsSy2 OIMSs represent a 
transition from SMGs to QSOs, and eventually to cur- 
rent epoch massive elliptical galaxies, it is a relatively 
rapid one. 

4.6. Do Interactions Trigger the Activity in 
OIMSs/ OIRSs? 

Evidence that the near environment of galaxies plays a 
role in triggering active nuclei has been known for some 
time. Hcckman ct al. (1984, 1986), for example, found 
that a large fraction of low redshift quasars have com- 
panion galaxies within ss 50 kpc that typically differ in 
radial velocity from the quasar by < 1000 km s" 1 . In 
addition, roughly one-third of radio galaxies show highly 
peculiar morphologies (i.e., tails, bridges, and shells) in- 
dicative of tidal interactions. More recent work supports 
these results (e.g., Canalizo & Stockton 2001). The 3.6 
/zm images in Figure 1 suggest that most IRAC detected 
OIRSs possess close companions, and may therefore be 
tidally interacting. If so, this implies that like low red- 
shift radio galaxies, the central activity in the majority 
of the IRAC detected OIRSs may be induced through 
collisions. 

The fact that most OIMSs appear point-like with 
IRAC might at first suggest that they are relatively 
isolated. However, the mid-IR images shown in Fig- 
ure 3 cannot exclude the possibility that OIMSs pos- 
sess faint or lower-mass companions capable of inducing 
tidal perturbations large enough to fuel an active nu- 
cleus. Nor can the IRAC images rule out the possibility 
that OIMSs are a population of late-state mergers. Local 
examples, such as NGC 7252, typically show luminous 
asymmetric cores with extended low surface brightness 
tails (Schweizer 1982). Neither would be detectable with 
IRAC given the OIMS's redshifts. It thus remains an in- 
teresting possibility that interactions have triggered the 
nuclear activity in both OIRSs and OIMSs. Deep and 
high angular resolution near-infrared imaging studies will 
be needed to test this hypothesis. 

5. CONCLUSIONS 

OIMSs can be split into two populations based on their 
rest-frame 1-10 (im SEDs and IRAC colors in concor- 
dance with their IRS spectra: most (12/16) are domi- 
nated by a heavily obscured radio-quiet active nucleus, 
while the remainder are powered by either a starburst or 
a composite starburst/active nucleus. The AGN dom- 
inated obsSyl and obsSy2 OIMSs in particular are ex- 
tremely luminous in the rest-frame near and mid-IR, with 
iaL„(5 fim) comparable to the most luminous local NL- 
RGs and quasars. They can be regarded as "buried" 
QSOs and likely represent the predecesors of current 

6 For example, none of the OIMSs in Houck05 and in the larger 
sample of OIMSs were detected by MIPS at 70 /im or 160 fim. 
Conversely, it had been previously noted that SMGs arc fainter 
than expected at 24 fim given their inferred Lir (e.g., Ivison et 
al. 2004; Egami et al. 2004). This difference may simply reflect 
the presense of a strong contribution from AGN heated dust in the 
OIMSs that is lacking in the "cold" SMGs. 



epoch massive elliptical galaxies. OIMSs are also dis- 
tinct from other high-z source populations routinely se- 
lected using UV/optical or optical/near- infrared criteria 
such as LBGs and BzKs in their dominant power source, 
levels of obscuration, and mid-IR luminosity. Compared 
with other optically faint mid-IR selected populations 
(e.g., Yan et al. 2007), OIMSs represent extremes in 
both obscuration and AGN luminosity. This follows from 
their respecive selection criteria, which for OIMSs favor 
steeper mid-IR continuua and higher levels of extinc- 
tion. There appears to be significant overlap between 
OIMSs and SMGs, with the sb/comb OIMSs appear- 
ing largely indistinguishable from SMGs in their optical, 
mid-infrared, and radio properties. Moreover, the ob- 
sSyl and obsSy2 OIMSs may represent a brief obscured 
phase in the transition of a "cold" SMG to a QSO, and 
eventually to a massive current epoch elliptical galaxy. 

The OIRSs do not represent a single source popula- 
tion. The minority (6/35) that are detected by both 
IRAC and MIPS at 24 /zm have SEDs, mid-IR colors, and 
q- values indicative of either starburst or composite star- 
burst/ AGN powered systems. For these sources we find 
•Zphot in between 1.0 and 4.5, implying rest-frame K-band 
luminosities, maximum star formation rates, and stellar 
masses virtually identical to the "cold" (i.e., starburst 
dominated) SMGs in the Lockman Hole East region. 

The remaining OIRSs - which comprise 83% of the 
parent population in Higdon05 - fall into two classes de- 
pending upon whether or not they are detected by IRAC. 
Those that are detected have flat mid-IR SEDs implying 
comparable luminosities from stellar photospheres and 
hot AGN illuminated dust. As a group, their average 
L24/ tm /L3.6^m ratio is most consistent with either Low 
Excitation Radio Galaxies (LERGs) or objects like the 
"red and dead" radio galaxy LBDS 53W091. Those that 
are not detected by IRAC must be at z > 2 if they are 
as massive as the hosts of local radio sources. Both of 
these groups are characterized by negative values of q 
and thus may represent a population of relatively unob- 
scured radio galaxies at high redshift. Both represent 
populations that are highly distinct from the obsSyl and 
obsSy2 OIMSs. 

Differences between the optically invisible populations 
detected through observations at sub-millimeter, mid- 
IR and radio wavelengths can be understood in terms 
of selection effects: (1) extreme optical/mid- IR lumi- 
nosity ratios appears to ensure highly obscured AGN 
dominated sources (obsSyl and obsSy2 OIMSs), (2) ex- 
treme optical/sub-millimctcr luminosity ratios will se- 
lect highly obscured sources primarily powered by star 
formation (SCUBA/MAMBO sources), and (3) compact 
sub-mjy radio sources lacking optical counterparts ap- 
pears to preferentially choose distant BL Lac-like objects 
(OIRSs). 

We find no evidence that OIMSs or OIRSs inhabit the 
cores of rich clusters. Nor do we find significant differ- 
ences in local galaxy density between the two on < 100 
kpc scales. However, unlike the OIMSs, a large fraction 
of the IRAC detected OIRSs appear to possess close and 
massive companions, though higher angular resolution 
studies will be needed to reach firm conclusions. This 
suggests that the luminosity of OIRSs (and conceivably 
OIMSs) may be triggered by tidal interactions, as ap- 
pears to be the case for low-z radio galaxies and quasars. 
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The IRAC detected OIRSs may thus represent the for- 
mation of very massive galaxies at high redshift through 
major mergers. 
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TABLE 1 

Observed Flux Densities For The IRAC Detected OIRSs 



OIRSM F Bw < b > 
[HHW2005] # (/xJy) 


(**Jy) 


Fj (b) 

(pJy) 


(**Jy) 


F4.5Mm 

(MJy) 


F5.8 Mm 

(MJy) 


Fg.OMm 

(^Jy) 




F24pm 

(mjy) 






F20cm 

(mjy) 






q (c) 


19 


<0 


.062 


<o. 


256 


<o. 


404 


2b. 7 ± l.y 


27.6 


± 


2.5 


25.0 ± 12.5 


<35.0 




<0.18 







.57 


± 


0. 


04 




<-0.50 


49 


<0 


.056 


<0. 


233 


<0. 


255 


12.5 zb 2.1 


15.6 


± 


2.7 


<40.5 


66.2 zb lz.l 




<0.18 







.57 


± 


0. 


03 




<-0.50 


79 


<0 


081 


<0. 


256 


<0. 


.193 


29.1 zb 2.1 


35.2 


± 


2.7 


<45.3 


43.9 ± 12.2 




<0.18 







.95 


± 


0. 


04 




<-0.72 


97 


<0 


.056 


<0. 


147 


<0. 


.212 


19.7 zb 1.9 


24.1 


± 


2.5 


<40.7 


on o _l_ 1 o 1 
30.2 ± 12.1 





23 i 0. 


01 





.10 


± 


0. 


.03 


-0 


o A _l_ n no 
.24 ± U.08 


114 


<0. 


.081 


<0. 


194 


<0. 


255 


16.8 zb 1.8 


18.4 


± 


2.3 


31.3 ± 12.6 


31.6 ± 11.2 




<0.18 







,37 


± 





05 




<-0.31 


156 


<0 


.067 


<0. 


122 


<0. 


176 


22.0 ± 1.9 


22.8 


± 


2.5 


<39.9 


<35.8 




<0.18 







.80 


± 


0. 


06 




<-0.65 


176 


<0 


.067 


<0. 


194 


<0 


336 


Q Q o _l_ o n 
38. 2 ± 2.2 


74.2 


± 


3.1 


i err q 1 t r -i 

155.3 ± 15.1 


O A O 1 _1_ 1 Q O 

245.1 ± 13.2 





.46 i 0. 


04 





.25 


± 





01 


0. 


2b ± 0.U8 


182 


<0. 


.067 


<0. 


161 


<0. 


.336 


y. / ± L.i 


8.7 


± 


1.8 


<oy.y 


<OO.U 




<U.lo 




u. 


.55 


± 





n 1 
UT 




^ n Act 

<-u.4y 


185 


<0 


.067 


<0. 


122 


<0 


212 


16.0 ± 1.9 


17.0 


± 


2.4 


31.0 ± 13.0 


<35.8 




<0.18 







33 


± 


0. 


04 




<-0.26 


208 


<0 


.067 


<0. 


122 


<0 


193 


32.1 ± 2.1 


35.8 


± 


1.7 


<40.3 


39.8 ± 12.1 





26 ± 0. 


08 





.19 


± 





03 


0, 


13 ± 0.15 


232 


<0 


.074 


<0. 


281 


<0. 


336 


8.1 ± 1.7 


9.1 


± 


2.3 


<40.4 


<35.9 




<0.18 




0. 


.51 


± 





02 




<-0.45 


245 


<0 


.074 


<0. 


256 


<0 


.306 


21.2 ± 1.9 


25.7 


± 


2.5 


41.5 ± 12.8 


25.0 ± 11.2 





44 ± 0. 


05 





,24 


± 


0. 


03 


0, 


26 ± 0.07 


278 


<0 


.076 


<0. 


148 


<0. 


.146 


<2.2 


6.1 


± 


1.6 


29.4 ± 13.5 


<35.6 




<0.18 







.15 


± 


0. 


07 




<-0.40 


346 


<0 


.067 


<0. 


112 


<0. 


.176 


8.8 ± 1.6 


13.2 


± 


2.2 


<36.8 


47.9 ± 11.4 




<0.18 







37 


± 


0. 


04 




<-0.31 


349 


<0 


.062 


<0. 


134 


<0 


232 


26.7 ± 0.9 


28.9 


± 


1.5 


22.5 ± 6.9 


38.3 ± 7.2 





.25 ± 0. 


08 





,24 


± 





05 


0. 


02 ± 0.17 


363 


<0 


.056 


<0. 


134 


<0. 


.193 


8.6 ± 0.9 


12.0 


± 


1.4 


<40.3 


<36.0 





.38 ± 0. 


05 





,15 


± 


0. 


04 





40 ± 0.13 


380 


<0 


.074 


<0. 


194 


<0. 


.279 


9.4 ± 1.7 


13.5 


± 


2.4 


<40.0 


<35.8 




<0.18 




1 


92 


± 


0. 


04 




<-1.0 


389 


<0 


.074 


<0. 


112 


<0 


193 


12.4 ± 1.8 


22.3 


± 


2.5 


<40.0 


<36.0 




<0.18 







30 


± 





02 




<-0.22 


410 


<0 


.067 


<0. 


147 


<0. 


.193 


28.2 ± 2.0 


32.5 


± 


2.6 


<40.3 


<36.0 




<0.18 







.88 


± 


0. 


.08 




<-0.69 



( a ) OIRS number from Table 2 in Higdon05, where source positions 
are given. These sources are also designated [HHW2005] by SIMBAD 
as (e.g., [HHW2005] 97).< b > 3<r limits from NDWFS imagcs>> q = 
log(F24 fim/F20 cm)- The radio flux densities are from Higdon05, while 
the 24 fim flux densities and limits arc from this work. 



TABLE 2 

Observed Flux Densities For The OIMSs 



[HSW2005] # (/iJy) 


(^Jy) 


Fj (b) 

(pJy) 


F3.6/jm 

(^Jy) 


F4.5 M m 

(MJy) 


F5.8(im 

(^Jy) 


F8.0(jm 

(^Jy) 




F24/jm 

(mjy) 




F 20 cm (c) 

(mjy) 


q (d) 


1 


0. 


169 


0. 


337 


0. 


443 


26.9 


± 


1.0 


46.0 


± 


1.4 


79.3 ± 


8.9 


207.7 ± 


8.2 


1 


.21 


± 





.01 


<0.15 


>0.90 


2 


<0 


.081 


<o. 


.102 


0. 


404 


17.0 


± 


1.0 


26.0 


± 


1.4 


20.4 ± 


8.0 


62.1 ± 


7.3 





89 


± 


0. 


01 


<0.15 


>0.76 


3 


<0 


.067 


<0. 


123 


<0 


.193 


22.3 


± 


1.0 


26.0 


± 


1.6 


39.0 ± 


7.4 


73.9 ± 


8.2 


1 


81 


± 





02 


0.31 ± 0.03 


0.77 ± 0.04 


4 


0, 


089 


0. 


445 


0. 


584 


17.6 


± 


0.9 


38.0 


± 


1.4 


92.4 ± 


6.4 


248.5 ± 


7.3 


1 


.08 


± 





01 


0.24 ± 0.04 


0.65 ± 0.06 


5 


<0 


.067 


<0 


.123 


<0 


.193 


4.8 


± 


0.8 


9.4 


± 


1.3 


<33.0 


<21.0 







87 


± 


0. 


01 


<0.15 


>0.75 


6 


<0 


.067 


<0. 


.123 


<0 


.193 


15.6 


± 


0.8 


21.1 


± 


1.5 


30.4 ± 


6.3 


92.6 ± 


8.2 


1 


03 


± 


0. 


02 


<0.15 


>0.82 


7 


0, 


141 


<0 


.123 


0. 


443 


41.7 


± 


0.8 


55.3 


± 


1.4 


70.1 ± 


6.5 


124.7 ± 


7.8 





78 


± 


0. 


01 


<0.15 


>0.70 


8 


0. 


129 


0. 


123 


0. 


279 


12.2 


± 


0.8 


18.8 


± 


1.4 


69.6 ± 


6.7 


231.3 ± 


7.9 


2, 


65 


± 


0, 


02 


<0.15 


>1.23 


9 


<0 


.067 


<0. 


.123 


<0 


.193 


7.4 


± 


1.0 


28.3 


± 


1.2 


92.4 ± 


6.8 


533.3 ± 


5.9 


3. 


83 


± 





02 


0.42 ± 0.03 


0.96 ± 0.04 


10 


0. 


511 


0. 


123 


0. 


336 


6.0 


± 


0.8 


12.5 


± 


1.4 


25.3 ± 


7.7 


<21.0 




1 


45 


± 


0. 


01 


<0.15 


>0.97 


12 


<0 


.067 


<0. 


123 


<0 


.193 


14.3 


± 


1.0 


25.9 


± 


1.3 


<41.0 


38.0 ± 


6.8 


1. 


.12 


± 





01 


0.20 ± 0.03 


0.80 ± 0.07 


13 


0, 


561 


0. 


929 


1. 


112 


42.1 


± 


3.2 


107.2 


± 


4.5 


328.1 ± : 


22.5 


699.2 ± 25.1 


2. 


30 


± 





01 


5.09 ± 0.03 


-0.34 ± 0.01 


14 


<0 


.067 


<0. 


123 


<0 


.193 


16.9 


± 


0.8 


22.6 


± 


1.5 


26.4 ± 


7.0 


25.7 ± 


7.8 


0, 


79 


± 





01 


<0.15 


>0.69 


15 


<0 


.067 


<0. 


123 


<0 


.193 


33.3 


± 


0.9 


48.5 


± 


1.4 


83.6 ± 


6.9 


148.8 ± 


7.5 


1. 


05 


± 





01 


<0.15 


>0.82 


16 


<0 


.067 


<0. 


.123 


<0 


.193 


3.6 


± 


0.9 


9.0 


± 


1.4 


16.0 ± 


7.5 


60.4 ± 


8.0 


1. 


.04 


± 


0, 


01 


<0.15 


>0.81 


17 


<0 


.067 


<0 


.123 


<0 


.193 


20.6 


± 


0.8 









56.2 ± 


6.3 


... (e 


) 


1. 


23 


± 


0, 


01 


<0.15 


>0.89 



OIMS number from Tabic 1 in Houck05, where source positions 
are given. These sources are also designated [HSW2005] by SIMBAD 
(e.g., [HSW2005] 13). ( b ) Optical flux densities or 3ct limits from 



NDWFS. (c) 20 cm flux densities or 5cr limits from de Vries et 
al. (2002). ( d > q = log(F 24Mm /F 20(; m). The -Fb^m values arc from 
Houck05.^ No data at these bands. 



14 



Higdon et al. 



TABLE 3 

Derived OIRS z phot And Luminosities 



OIRSM 


Zphot 


Template' 13 ) 


AyM 


L K < d > 


1 r ( e ) 




L20cm 


[HHW2005] # 








( 10 11 L 0iK ) 


( 10 37 W ) 


( 10 11 M ) 


( 10 25 W Hz" 1 ) 


97 


1.5 ± 0.3 


Arp 220 


2 


3.7 ± 1.6 


0.3 ± 0.2 


1.4 ± 0.5 


0.4 ± 0.2 


176 


4.5 ± 0.2 


M 82 


1 


278. ± 26. 


11.6 ± 1.3 


116. ± 11. 


2.9 ± 0.3 


208 


2.0 ± 0.3 


M 82 


1 


9.9 ± 3.4 


0.7 ± 0.2 


3.4 ± 1.1 


0.4 ± 0.1 


245 


2.0 ± 0.3 


Arp 220 


2 


7.4 ± 2.7 


1.2 ± 0.3 


3.3 ± 1.1 


0.5 ± 0.2 


349 


1.0 ± 0.2 


M 82 


2 


1.8 ± 0.6 


0.09 ± 0.03 


0.8 ± 0.3 


0.2 ± 0.1 


363 


2.2 ± 0.6 


Arp 220 


1 


4.5 ± 2.4 


1.1 ± 0.6 


2.3 ± 1.4 


0.4 ± 0.2 



( a ) OIRS number from Table 2 in Higdon05, where source positions 
are given. These objects are designated in SIMBAD as [HHW2005] 
# (e.g., [HHW2005] 97). ( b ) The best fitting galaxy template; see 
§3.1 for descriptions. ( c ) Additional visual extinction in magnitudes 
to satisfy optical constraints.^) The source rest-frame Johnson Ri- 
band luminosity, defined as Lk — 47rD 2 (1+z) - 1 .FV, Ai^k, in units of 
the sun's K-band luminosity, L 0jK = 8.7 x 10 24 W (Bcsscl 1979). D L 
is the luminosity distance, Az/k is the Johnson K-band filter FWHM 
(2.1 X 10 13 Hz; Tokunaga et al. 2002), and F v is the source flux 
density at A obs = (1 + z)2.2fim estimated from the source's SED 
template.^) The source rest-frame monochromatic luminosity at 5 fim 
in Watts, defined as i/L, y (5 jum), where L^(5 /im) — 47tDl (1+z) ~~ 1 F u . 
Dl is again the luminosity distance and F v is the source flux density at 
Aobs — (l-\-z)5fj,m estimated from the source's SED template.^) Stellar 
mass derived by multiplying Lk by M*/Lk — 0.3, which is appropriate 
for a starburst galaxy (Bell &c dc Jong 2001). 



TABLE 4 
Derived OIMS Luminosities 



OIMSM 


2spec 


Template^) 




L 5Mm < d ) 


M„( c ) 


L20cm 


log L IR ( f ) 


[HSW2005] # 






( 10 12 L , K ) 


( 10 38 W ) 


( 10 11 M ) 


( 10 25 W Hz" 1 ) 


(L ) 


1 


2.64 ± 0.25 


obsSy2 


8.3 ± 1.8 


8.3 ± 1.8 


<25. 


<0.8 


13.3 


2 


1.86 ± 0.07 


sb 


0.7 ± 0.1 


1.0 ± 0.1 


1.9 ± 0.2 


<0.4 


13.0 


3 


1.78 ± 0.30 


obsSyl 


0.9 ± 0.3 


0.9 ± 0.3 


<2.7 


0.5 ± 0.2 


13.0 


4 


2.59 ± 0.34 


obsSy2 


9.4 ± 2.9 


9.4 ± 2.4 


<28. 


0.9 ± 0.3 


13.2 


5 


2.34 ± 0.28 


obsSy2 








<0.7 


13.0 


6 


1.96 ± 0.34 


obsSyl 


1.2 ± 0.5 


1.7 ± 0.6 


<3.7 


<0.4 


12.8 


7 


1.78 ± 0.14 


comb 


1.6 ± 0.3 


1.4 ± 0.3 


4.8 ± 0.9 


<0.3 


13.4 


8 


2.62 ± 0.26 


obsSy2 


9.1 ± 2.1 


10.2 ± 2.2 


<27. 


<0.8 


13.6 


9 


2.46 ± 0.20 


obsSy2 


16.3 ± 3.2 


21.0 ± 4.1 


<49. 


1.4 ± 0.3 


13.7 


10 


2.08 ± 0.21 


obsSyl 


0.9 ± 0.2 


0.8 ± 0.2 


<2.8 


<0.5 


13.0 


12 


2.13 ± 0.09 


comb 


0.7 ± 0.1 


0.4 ± 0.1 


2.1 ± 0.2 


0.5 ± 0.1 


13.8 


13 


1.95 ± 0.17 


obsSy2 


10.2 ± 2.1 


10.9 ± 2.2 


<31. 


10.2 ± 2.1 


13.3 


14 


2.26 ± 0.11 


comb 


0.8 ± 0.8 


3.0 ± 0.4 


2.4 ± 0.3 


<0.6 


13.7 


15 


1.75 ± 0.21 


obsSyl 


1.8 ± 0.5 


1.6 ± 0.4 


<5. 


<0.3 


12.7 


16 


2.73 ± 0.19 


obsSy2 


2.7 ± 0.4 


2.9 ± 0.5 


<8. 


<0.8 


13.2 


17 


2.13 ± 0.17 


obsSyl 








<0.5 


13.0 


( a ) Source number from Table 


1 of Houck05 


where coordinates 


can 









be found. The OIMSs can also be referenced by their SIMBAD 
designation HSW2005 # (e.g., [HSW2005] 13) (b) Best fitting template 
spectrum type; see §3.2 for descriptions.^ Rest-frame Johnson K-band 
luminosity in units of the sun's K-band luminosity, Lq ( k — 8.7 x 10 24 
W (Bessel 1979). ^ Monochromatic luminosity at rest-frame 5 fim in 
Watts, defined as v\^ v (5 ,um).( e ) Stellar mass derived by multiplying 
Lk by M*/Lk — 0.3, which is appropriate for a starburst (Bell & de 
Jong 2001). ( f ) Estimated 8-1000 ^m luminosity from Houck05. 
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TABLE 5 

Sub-millimeter Galaxies In The Lockman Hole East Region 



Source w 


Z W 


TypeW 


L K (d) 


L 5M m (G) 




L20cm 


qis) 


[LE850] # 






( 10 11 L ,k ) 


( 10 37 W ) 


( 10 11 M Q ) 


( 10 24 W Hz" 1 ) 




1 


2.6 


cold 


4.0 ± 1.0 


1.4 ± 0.3 


1.2 ± 0.3 


3.1 ± 0.7 


0.42 ± 0.10 


4 


2.6 


cold 


3.8 ± 1.0 


1.3 ± 0.3 


1.1 ± 0.3 


1.3 ± 0.3 


0.68 ± 0.51 


7 


1.8 


cold 


12.2 ± 4.7 


4.2 ± 1.5 


3.7 ± 1.4 


2.0 ± 0.8 


0.36 ± 0.07 


8a 


1.0 


cold 


5.0 ± 3.5 


1.8 ± 1.0 


1.5 ± 1.1 


1.2 ± 0.7 


1.73 ± 0.24 


8b 


3.0 


warm 


53.2 ± 12.5 


18.4 ± 4.3 


<16. 


4.1 ± 1.0 


0.69 ± 0.14 


8c 


0.9 


cold 


1.4 ± 0.6 


0.6 ± 0.3 


0.4 ± 0.2 


<0.3 


>0.81 


14a 


2.4 


cold 


6.7 ± 2.0 


2.3 ± 0.7 


2.0 ± 0.6 


0.7 ± 0.3 


0.36 ± 0.24 


14b 


2.5 


cold 


6.9 ± 1.8 


2.4 ± 0.7 


2.1 ± 0.5 


1.6 ± 0.6 


0.36 ± 0.15 


18 


2.7 


warm 


11.4 ± 3.0 


3.9 ± 1.0 


<3.4 


1.1 ± 0.3 


0.42 ± 0.16 


35 


3.0 


cold 


20.3 ± 4.8 


7.0 ± 1.7 


6.1 ± 1.4 


2.8 ± 0.7 


0.44 ± 0.13 



( a ) Source number in Tabic 1 of Egami ct al. (2004), where radio, 
mid-IR, and sub-millimeter flux densities are given. Coordinates 
can be found via SIMBAD using the designation LE850 followed 
by the source number (i.e., [LE850] 4).( b ) Sources 8a, 14a, and 18 
have spectroscopic redshifts. The other redshifts are photometric. 
See footnotes a and b in Tabic 1 of Egami et al. (2004) for 
references/ ) Sources classified as "warm" possess power-law SEDs 
in the Spitzer data and are likely AGN powered, whereas "cold" 
sources appear to be powered largely by star formation (e.g., 
bottom-right panel of Figure 7).( d ' Rest-frame Johnson K-band 
luminosity in units of the sun's K-band luminosity, Lq k = 8.7 X 

10 24 W (Bessel 1979). (°) Monochromatic luminosity at rest-frame 
5 (jrn in Watts, defined as ^L„(5 (im).' f ' Estimated stellar mass 
from multiplying Lk by M*/Lk = 0.3, which is appropriate for a 
starburst (Bell & de Jong 2001). (s) q = ^(i^m/i^Ocm). 
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Fig. 1. — Optically Invisible Radio Sources (OIRSs) from 
Higdon05 with detections in the IRAC Shallow Survey (Eisenhardt 
et al. 2004). Optical Bw images are shown in the left-most column 
along with source index and estimated photometric redshift (if 
known) as discussed in §3.1. The middle-column shows the same 
field with IRAC at 3.6 /im, along with 20cm radio contours starting 
at 3 a from Higdon05. The third column shows the IRAC 3.6/im 
image without radio contours. The fields are 20" wide, which 
corresponds to «175 kpc for z = 1 — 3. The source index refers 
to Tabic 2 of Higdon05 and Table 1 of this paper. A linear stretch 
is used for both optical and infrared grey-scale images. 



Fig. 2. — Optically Invisible Radio Sources (OIRSs) with no 
obvious counterparts in the IRAC Shallow Survey (Eisenhardt ct 
al. 2004). Optical Bw images are shown in the left-most column 
while the corresponding IRAC 3.6 /im images are shown to the 
immediate right. Radio 20cm contours are shown in the middle 
IRAC 3.6 fim planes, starting at the 3 a level. The fields are 20" 
in size. Also shown are the source index from Table 2 of Higdon05 
and Table 1 of this paper. A linear stretch is used for both optical 
and infrared grey-scale images. 



Fig. 3. — Optically Invisible 24/im Sources (OIMSs) from 
Houck05. Optical Bw images arc shown in the left-most column 
along with rcdshifts derived from IRS spectra. Opposite these are 
matching IRAC 3.6/im images from the IRAC Shallow Survey 
(Eisenhardt ct al. 2004). Each field is 20" wide, which corresponds 
to Ril75 kpc for 2 = 1 — 3. The source index refers to Table 1 of 
Houck05 and Table 2 of this paper. A linear stretch is used for both 
optical and infrared grey-scale images. 
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Fig. 4. — Minimum-x 2 template fits for the comb and sb 
OIMSs. The Arp 220 SED template (solid line) gives photometric 
rcdshifts in good agreement with IRS derived z Bpec (Az « 0.15). 
Wavelengths are shown in the observer's frame on the bottom axis 
and in the rest-frame on the top axis. 
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Fig. 5. — Minimum-x 2 template fits for the six IRAC and MIPS 
detected OIRSs. The galaxy template that best fits the data and 
the derived Zphot are listed at the top of each panel. Upper-limits 
arc 3<t, and arc depicted as arrows. Wavelengths are shown in the 
observer's frame on the bottom axis and in the rest-frame on the 
top axis. 
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Fig. 6. — Minimum-x 2 template fits for two OIRSs detected 
to > 3cr in at least three IRAC bands but undetected by MIPS. 
The axes are defined as in Figures 4 and 5. The derived z p hot m 
each instance is only an estimate. Similarly reasonable fits using 
an elliptical template are depicted using a dashed line. For either 
template choice, the estimated z p hot does not change substantially. 
Wavelengths are shown in the observer's frame on the bottom axis 
and in the rest-frame on the top axis. 
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Fig. 7. — Averaged rest-frame SEDs for the OIMSs grouped 
by the template galaxies that best fit their IRS spectra: obsSy2 
(IRAS F00183-7111), obsSyl (Mrk 231), comb (Arp 220), and sb 
(NGC 7714). Also shown are averaged rest-frame SEDs for the 
six IRAC detected OIRSs with z p ^ ot and the eight "cold" SMGs 
from Egami et al. (2004) and Ivison ct al. (2004). Each SED has 
been normalized to A ros t = 1.4 fim. Individual data points (limits) 
are represented by dots (small arrows) while the unfilled squares 
represent binned averages. Large arrows represent averaged limits. 
Note that the sft-likc SED in the middle-right panel consists of a 
single source (OIMS #2) and has not been rebinned. 
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Fig. 8. — IRAC two-color diagrams. The top panel shows 
OIMSs (filled squares if their IRS spectra are obsSyl or obsSyS- 
like and empty squares if sb or comfe-like) along with ten SMGs 
in the Lockman Hole East region (asterisks for "cold" and six- 
pointed stars for "warm" sources) from Egami et al. (2004). Tracks 
corresponding to Mrk 231 and a dusty M 82 template (extra Ay 
= 2.5) for 2 = 1 — 5 are over-plotted. OIMS #17 was not observed 
at 4.5 fim and 8.0 fim and is not shown. The lower panel shows 
the OIRSs from Table 1 as circles. Those with z p hot are unfilled. 
A track corresponding to a model elliptical galaxy for z = 1 — 5 is 
also shown. 
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Fig. 9. — The rest-frame -F24nm/F3.6/jm ratio as a function of 
redshift for five sources known to contain powerful AGN: (1) the 
averaged obsSy2 SED shown in Figure 7, (2) the NLRG 3C459, (3) 
the quasar 3C273, (4) the LERG 3C293, and (5) the "red and dead" 
radio galaxy LBDS 53W091 (Stern et al. 2006). The solid horizontal 
line at .F24|tim/-F3.6j«n = 3.9 indicates the 3<r upper-limit for the 
twelve IRAC detected OIRSs without 24 fim detections, defined by 
dividing the stack-averaged upper-limit at 24 fim by the average 
3.6 fim flux density. The two horizontal lines above and below 
this show the range in F , 24/jm/P3.6/jm for these OIRSs defined by 
dividing the 24 fim upper-limit by their individual 3.6 fim flux 
densities. These OIRSs appear to have SEDs most consistent with 
LERGs and "red and dead" radio galaxies like LBDS 53W091 over 
a wide redshift range. 
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Fig. 10. — A radio-infrared two-color diagram showing OIMSs, 
IRAC detected OIRSs, and the ten SMGs in the Lockman Hole 
east region observed by Egami et al. (2004). The AGN-dominatcd 
OIMSs (i.e., obsSyl and obsSy2) are depicted as filled squares 
while the OIMSs with a significant starburst component {comb 
and sb) are depicted as empty squares. The OIRSs are shown 
as unfilled circles if they have z p hot estimates, and filled circles 
if no redshift estimates are available. Similarly, "cold" starburst 
dominated SMGs are shown as asterisks, while the "warm" AGN 
dominated SMGs are shown as six-pointed stars. The dashed 
horizontal lines show the limits of the radio-infrared correlation 
determined by Appleton et al. (2004) at 24 fim and 20cm (q = 0.8 
± 0.3). OIRS #278 and OIMS #17 could not be plotted. 
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Fig. 11. — Rest-frame K-band luminosity as a function of 
observer-frame Fg. o^m/ Fss/im flux density ratio for OIMSs (filled 
squares for obsSyl and obsSyS-likc systems, empty squares for sb 
and comfe-like systems) and ten SMGs from the Lockman Hole 
East region observed by Egami et al. (2004). Starburst-like "cold" 
sources are shown as asterisks while AGN-like "warm" sources are 
depicted as six-pointed stars. The six IRAC detected OIRSs with 
photometric redshifts are also shown as empty circles. 
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